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Part I. 

§ 1. On the Scattering of Parallel Radiation by Molecules and Small Particles. 

The effect of small particles in scattering incident radiation was first worked out by 
Lord Rayleigh.^ When a stream of parallel radiation falls on a particle whose 
dimensions are small compared with the wave-length the resulting secondary 
disturbance travels in all directions at the expense of the intensity in the original 
direction. In a later paper Lord RAYLEiOHf gave reasons for believing that the 
molecules of a gas are themselves able to scatter radiation in this way. In a gaseous 
medium it is legitimate to sum up the intensities of the scattered radiation due to 
each molecule in an element of volume without a consideration of phase-difference in 
consequence of the continuous change in the relative positions of a molecule in a gas. 

^ Eayleigh, 'Phil. Mag.,' XLL, pp. 107, 274, 447 (1871); XII., p. 81 (1881); 'Collected Works,' 
vol. I, pp. 87, 104, 518. 

t Eayleigh, 'Phil. Mag.,' XLVIL, pp. 375-384 (1899); ' Collected Works,^ vo}. IV., pp. 397-405. 

VOL. GGXII. A 495. Publislied separately, February 8, 1913. 
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The same remark applies to the case where the scattering is due to small particles of 
dust since these partake, to some extent at least, of the molecular agitation of the gas 
in which they are held in suspension. 

The intensity of- radiation at any point is measured by the amount of energy 
crossing unit area of a surface normal to the direction of the radiation in unit time. 
Unpolarized radiation of intensity E falls on an element §v of a gas of density p and 
containing N molecules per unit volume. The intensity of the scattered radiation in 
a direction 6 with the incident radiation and at a distance r from §v we denote by 
r^^l {r, 6) Sv, so that I (r, 6) §v Sw is the energy contained in a small solid angle §m 
crossing a spherical surface at distance r in unit time. 

The expression for I (O, 0) is of the form 

I(O,0) = m{6)E . . . . . . . . . . (1) 

where jul (d) depends on the direction 6 and is proportional to the number N of 
molecules per unit volume, i,e., /x{0) is also proportional to the density p. 

If /Ao (O), No, pq refer to values of fi (6), N, p under determinate conditions of pressure 
and temperature we have 

M (0)/mo (e) - N/No - p/po (2) 

JUL (d) may be expressed in terms of the optical properties of the medium : the results 
of Rayleigh,^ and KsLViNt worked out on various hypotheses of the molecule and of 
the aether agree in giving rise to the expression, 

^(d) = i7r'{n'-iyX-'{l-^C0B'9)/^ ....... (3) 

where n is the refractive index of the gas and X the wave-length of the incident 
radiation. 

Schuster| has recently obtained this result from general considerations indepen- 
dently of any particular theory. 

Since n—l and N are both proportional to the density of the gas we notice that 
JUL (6) is also proportional to the density, as already stated. 

Formula (3) shows that the intensity of the scattered radiation is twice as great 
in the direction of the incident radiation as it is in a direction at right angles. 
Considerable simplification in the analysis is obtained by employing the mean value of 
^a{d) over a spherical surface. Thus, writing 

4i7rfl = j iJL\u)dw .......... (4) 

(3)^gives 

fi = %^^{n^-^iy\-'l^ ......... (5) 

■^ Eayleigh, *Phil. Mag.,' XLI., p. 107 ; ' Collected Works,' vol. I., p. 87. 

t Kelvin, * Baltimore Lectures' (1904), p. 311. 

X Schuster, ^ Theory of Optics/ 2nd ed. (1909), p. 325. 
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and under standard conditions of pressure and temperature, 

flo = W{no'-iyX-i^, (6) 

We denote by kq the expression 

K, = iTPi, = %7rHn,'-iy\-i'^, (7) 

The effect of scattering is to diminish the intensitj^ of the incident radiation, giving 
rise to the phenomenon of attenuation especially noticeable in the diminution of 
intensity of solar radiation in its passage through the earth's atmosphere. The 
consideration of attenuation as due to scattering alone involves the assumption that 
energy is nowhere accumulating in the gas. In order to give greater generality to 
the application of the analysis we introduce a term expressing the fact that the 
temperature at any point is increasing. 

If E be the intensity of radiation crossing unit area of a plane at a point x in unit 
time, the loss to E in a distance dx in unit time due to the conversion of radiant 
energy into molecular agitation represented by a rise of temperature is of the form 

dE = -aEdx, • (8) 

when a is proportional to the number of molecules per unit volume, i.e., if ao and pq 
refer to standard conditions of pressure and temperature, a/oto = p/po. 

In the case of a pure gas a is a quantity depending on the distribution of energy 
between vibrating systems within the molecules and the motions of the molecules 
themselves which define the temperature of a gas on the kinetic theory scale as 
proportional to the mean squares of molecular velocities. 

If p be the density of the gas, s its specific heat at constant pressure and dO the 
increment of temperature in time dt due to the conversion of radiant energy aE dx dt 
into heat, we have 

aE = p5 ^j , (9) 

E being measured in calories per unit area normal to the direction of E per second. 

In the case of a gas exposed continuously to external radiation we may suppose the 
temperature to attain to a steady state when dO/dt = and therefore a = in which 
case no energy accumulates in the medium. In the problem of the earth's atmosphere, 
however, the term a will give rise to a small diurnal variation of temperature 
throughout the atmosphere. The term a may also be taken to include the effect of 
dust in absorbing solar radiation without scattering as well as effects of selective 
absorption if we regard a as a function of the wave-length. The existence and 
magnitude of this quantity can be determined by a comparison of theoretical results 
with the results of observation. Actual numerical values for air are obtained in 
Part III. of the present paper. 

VOL. ccxii. — A. 3 o 
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§ 2. On the General Integral Equation for the Scattered Radiation, 

Suppose a mass of gas to be completely enclosed by a boundary S and illuminated 
by a distribution of radiation whose intensity at a point {x, y, z) is E, and whose 
direction is defined at that point. Each element of volume will scatter a certain 
proportion of the radiation incident upon it, so that each element besides being 
illuminated by the incident radiation is also subject to the aggregate radiation from 
all the other elements within the surface 2, i.e., to the effect of self-illnmination. 
This constitutes what Schuster'^ has called the problem of '^ Radiation through a 
Foggy Atmosphere/^ The problem does not lend itself easily to a complete formula- 
tion in terms of differential equations, but can be expressed in terms of an integral 
equation. 

At an element of volume 8v at {x^ ?y, z) the incident radiation is E. At a distance 
r from 8v in a direction Q with the incident radiation the scattered radiation crossing 
unit area in unit time is denoted by I (r, Q) Sv r~^. If we wish to express the 
dependence of I (r, 6) Sv r~^ on the position of Sv in the surface 2 we write it in the 
form I {x^ y^ z^ r, Q) 8v r~^. 

Our first problem will be to express the scattered radiation I (r, Q) in terms of the 
scattered radiation I (O, Q) per unit solid angle in the neighbourhood ^v. 

Consider the radiation I (r, B) Sv Soo coming from an element of volume Sv and 
contained in a small solid angle Sm in a direction 6 with the incident radiation. The 
intensity of the radiation crossing a spherical surface of radius r is I (r, 6) Sv Sw ; that 
cutting the surface r + ^r is 



I (r, 0) + j^ I (r, 6) Sr \ Sv Sw. 



We thus obtain the equation 



- ^ I (r, e) Sv Sw = al {r, 6) r"' Sv r^ Soy c^r + 4^/iI (r, 6) r~' Sv r' Sw Sr, . . (lO) 

or 

The first term on the right-hand side represents by (8) the energy lost to the element 
of volume r^ Sw Sr by the conversion of radiant energy into a rise of temperature or 
into long- wave heat radiation, which transformation we assume to go on at a 
constant rate. 

The second term accounts for the energy lost to the element of volume r^ Sw Sr by 
scattering. We neglect the effect of self-illumination within the small solid angle Sai 
since this only enters into the equation to the second order of small quantities. 

■^ Schuster, "Radiation through a Foggy Atmosphere," ' Astrophysical Journal,' XXI, January, 1905, 
p. 1. A somewhat similar problem had been previously considered by the same writer in a paper ^' The 
Influence of Radiation on the Transmission of Heat," 'Phil. Mag.,' February, 1903. Recently it has been 
shown by Jackson, W. H. ('Bull. Am. Soc. Math.,' XYI, June, 1910, p. 473), that the generalized 
differential equation obtained in Schuster's paper can be transformed into an integral equation of the 
Fredholm type. 
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Writing k — iirjl the above equation gives 

-^J{r,0) = (a + /c)I(r,0) , . . (ll) 

If we write 

= OL-\-Kj \^^} 

the solution of (IX) is 

I(r,0) = I(O,0)6-^>'% . ........ (13) 

which puts into evidence the variation of a and k with the density of the gas. 
Consider now the radiation incident on an element of volume h at (cc, y, z) : — 

(i) The external illumination E which contributes to the scattered radiation 

from Sv an amount 

JUL {$) E {x, y^ z) Sv Sw ; 

(ii) An element of volume Sv^ at {x\ y\ z^) gives rise to an intensity 
I {x\ y\ z\ 7'\ 6^) Sv^ r^~^ at the point {x, y, z) : this contributes to the scattered 
radiation from Sv the amount 

M (P) I {x\ y\ z\ r^ 00 ^^' ^^ ^^ ^^'"'. 

TT^ denotinff the angfle between r and t\ The total contribution to the 
scattered radiation from Sv due to self- illumination by the entire volume 2 is 

So, Sv \ M (P) I (^', y\ z\ r^ 00 ^^'"' ^^'^ 

the integral being taken throughout the entire volume enclosed by the 
surface 2. 

By definition the sum of contributions (i) and (ii) is 1(^,^,0,0,0), so that we 
obtain the following integral equation for the scattered radiation at and from any 
point, 

l{x,y,z,Q,e) = ^ie)~E{x,y,z)+ f p,(^^)l{x\y\z\Q,ff)r'"'e-^>'''Uv\ . (14) 

I {x\ y\ z', r\ 00 being expressed in terms of I {x\ y\ z\ 0, 00 by means of (13). 

A differential equation involving E {x, y, z) is obtained by considering the rate of 
accumulation of energy in an element of volume. 

As soon as I {x, y, z, 0, 0) is known as a function of the position of the point O, 
{x, y, z), in the volume bounded by the surface 2, the radiation scattered to any point 
P within the boundary 2 contained in a small solid angle m is given by the formula 

Tw = w I(x,^, 0, 0, 0)6""J«^^'^''(ir, (15) 

Jo 

where r = PO, r^ - PQ, the radius vector to the boundary, and T is the intensity at 

Q Kj ^ 
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P for unit solid angle. The formula (15) is easily modified to include the case where 
P is without the surface 2. 

It is interesting to point out the analogy which the method of the present section 
bears to the ordinary procedure of potential theory. The function I {x, y, z) 
corresponds to a potential function and is expressed in terms of an external effect 
E {x, 2/, z) by means of an integral equation. The total intensity in a small solid 
angle, the effect which is physically measurable, is derived hj integrating between 
definite limits along a given direction. 

Further progress towards a solution of (14) is impossible without a number of 
simplifying assumptions which are best considered in dealing with a particular 
problem, such, for instance, as that presented by the effect of the earth's atmosphere 
in absorbing and scattering the sun's radiation. 



Part II. 

§ 3. Application to Radiation and Absorption in the Earth's Atmosphere, 

In the following sections we shall assume as an approximation that the surface of 
the earth is a plane and that the density is a function of the height above the earth's 
surface only. We also neglect effects due to reflection from the earth's surface and to 
refraction by the earth's atmosphere. 




'entth 



Sv* 




Fig. 1. 

The integral equation (14) for the scattered radiation in a direction (fig. l) with the 
direction of the sun and coming from an element ^v of the atmosphere at a height x 
above the earth's surface, takes the form 



l{x,Q,e)-= fx{e)'K{x)^\ fx{rr') 



qj2 



dv, . 



. . . (16) 



where the integral is taken throughout the entire atmosphere. 

In order to simplify the solution of (16) we assume that the scattered radiation 
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from a molecule can to a first approximation be taken equal in all directions, i.e., the 
equation can be reduced to one with a single variable by writing 

JUL (O) ~ JUL (rr^) — fl. 

Equation (16) may now be written, since dv' = r^'^ dod dr\ 

i(^)-mE(^)+/zJ 1(^0^"^"''''^^'^^^' • • W 

Eemembering that jul/juq = k/kq = a/^o = K/Ko = p/po, we now employ the transforma- 
tions 



.J, r.... ^CC 



n=\^dr, X=\^dx, H= £-dx, J{X) = ^I{x),, . . (18) 

Jo po Jo po Jo ^0 P 

and (17) now takes the form 

J(X) = MoE(X) + Mjj(X')e-'^"«'c?a,'c?E', (19) 

while (15) transforms into 

Tft, - w r° J (X) e-^»« dB (20) 

Jo 

The integral on the right-hand side of (19) must now be taken throughout a 
homogeneous atmosphere of density p^) included between the planes X = and X = H. 
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Fig. 2. 
Expressing the integral in (18) in cylindrical co-ordinates {f, \jr') (fig. 2), we have 

where \fr^ is the azimuth of the element of volume dV^ referred to some fixed direction 
of reference. The integral then becomes 



j J (XO c?X'|" dir' ^y^^^r'.ir-xn'i^_lM. 



i"+{x'-xy 

The last integral when integrated with respect to i' gives 



• CO 

e 





-K„u-+a'-xrfi' i'J i' - = -Ei { -Ko (X'-X)} when X<X'<H 

^'^+(X— X)^ 

and = -^^•{-Ko(X-X')} when 0<X'<X, 
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where Ei{—x) is Glaisher's"^ Exponential Integral denoted by 



/.OO 



Ei{—x)-=^— u'^^e'^^'du. (21) 



*J X 



The mtegral equation (19) now takes the form 



"rx 



il 



J(X) = MoE(X)-27rMo J{X')Ei{-K,{X-X')}dT+ J{X')M{-K,{X'~X)}dX' . (22) 

The differential equation by means of which E may be expressed as a function of 
X can be obtained by a consideration of the rate of accumulation of energy of the 
direct solar radiation between the planes x and x + dx. An analysis similar to that 
by means of which equation (13) was obtained, together with the transformations of 
equation (18), lead to the expression 

E (X) = Se^"^^"«>-^^^; .....,.., (23) 

where S is the intensity of solar radiation outside the earth's atmosphere corresponding 
to a given wave-length and ^ is the zenith distance of the sun. 

If we consider the rate of accumulation of the total energy (including both direct 
and scattered radiation) between the planes x and x + dx, it can be shown by making 
use of the integral equation connecting the direct with the scattered radiation that 
the exponential law of transmission expressed in (23) is valid, t] 

It is well to state clearly the assumptions involved in obtaining (13), (19), and (23) :— 

(i) In obtaining the differential equation leading to (13) the direct radiation is 
considered independently of the scattered radiation within the small solid 
angle w. 

(ii) The integral equation (22) assumes as an approximation that the radiation 
scattered by an element of volume is distributed equally in all directions. 

(iii) With these two conditions it can be shown from a consideration of attenuation 
of total radiation in a thin layer dx that the ordinary exponential law of 
transmission (23) follows ; i.e.^ that the transmission of direct radiation 
may be considered independently of the scattered radiation. 

(iv) By means of the transformation (18) it is shown that the problem of scattered 

* Glaisher, 'Phil. Trans.,' 1870, p. 367. 

[t JSfote added September ^0, 1912. — The calculation referred to was given at length in the paper as 
originally communicated ; the analysis is, however, somewhat tedious and hardly necessary in view of the 
fact that (13) is obtained according to the assumption that the direct radiation is considered independently 
of the scattered radiation within the small solid angle w. It is to be expected, therefore, that the 
attenuation of direct radiation in a parallel beam of solar radiation may be considered independently of 
the scattered radiation consistently with the above assumption leading immediately to (23) ; the considera- 
tion of the attenuation of the total radiation confirms this point and leads to the same result. The writer 
is indebted to the referees for the above suggestion.] 
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radiation in the earth's atmosphere reduces to the case of an atmosphere of 
uniform density contained between two parallel planes X = and X = H. 
This transformation is independent of any law of density with height, 
provided the planes of equal density are parallel to the earth's surface. 

§ 4. On the Approximate Solution of Integral Equations. 
The integral equation (21) is of the Fredholm type,"^ 

u{x) ■=f{oc)-{- ti{i)K{x,^)di, ....... (24) 

Except for special forms of the kernel K {x, i) the formal solutions of (24) do not 
lend themselves easily to numerical evaluation : we therefore develop a method of 
approximation which applies with sufficient accuracy to the problem in hand. 

Suppose for all values of x between Xi and x^ that /(a?) lies between A and a (A>a). 

Then to a first approximation 

u {x) lies between A + A K{x, i) di and a+ a K{x,^) d^, 



Xl 



Xi 



provided K (a^, ^) is everywhere positive. 
We write 

(f>{x)= K.{x, ^)d^, (25) 



X/i 



then if for all values of x, (p (x) lies between B and h (B>&) we have to a second 

approximation 

a + (a + afe) 5 <t^ (x) < A + (A + AB) B, 

or 

a{l-^b-\-¥)<u{x)<A{l-\-B + B'). 

A repetition of the process shows that 

If I B I < 1 both series are convergent and 

€i<u{x)<e2 where €i = a/l—b and e^ = A/l— B. . . . (26) 
Substituting in (24) we see that the solution u (x) lies between the limits 

Ui{x) =f{x)-\-€i(l){x), and U2{x) =f{x)-{'€2(p{x). .... (27) 
Uy (x) and U2 (x) may be called the extreme solutions. 

* Cf. BocHER, M., 'An Introduction to the Study of Integral Equations,' Cambridge, 1909, p. 14; 
also Bateman, H., ''Eeport on tlie History and Present State of the Theory of Integral Equations," 
' Brit* Assoc. Report,' 1910, p. 25. 
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If now a represent the mean value of /"(a;) between x^ and tCj, i.e. 



1 



a 



2 ] 



T yXr) Cl/ Jy J« • • • • • • • \-^^/ 



X\ 



and ^ stand for the mean value of ^ {pc) between x^ and x^^ i.e., 



/3 



1 






f/> {x) dx^ . 



x^ 



(29) 



while e = a/(l— /3), the solution u{x) =f{x)-\-€(p{x) may be called the mean solution, 
while the three solutions may be expressed in a single formula by the notation 



ic{x) =/(a?) + | € )^.{x) 



a « 



(30) 



We notice that Ui (x) < u (x) < U2 (x). In the applications to be considered, Ui (x) 
and U2 (x) are sufficiently close to warrant the use of the mean solution u {x). As 
long as u {x) does not depart far from the arithmetical mean of the extreme solutions, 
|- {ui (x) -{- U2 {x)} , we may take the value of u (x) to represent an approximation not 
far removed from the exact solution of the integral equation (24). 

It is perhaps worth noting here that to a higher degree of approximation the 
approximate solution of the integral equation may be written 



he^ 



u{x)==f{x)+ f{i)K{x,i)di+( ^e )i>{x),, .... (31) 

In the following sections the solutions (31) involve the evaluation of troublesome 
integrals, so that the simpler but somewhat less accurate approximations given in (30) 
will be employed. 



§ 5. On the Solution of the Integral Equation for Sky Radiation, 

In the integral equation (22) we may express to a first approximation the 
dependence of the scattered radiation in any direction on the angle which that 
direction makes with the incident radiation by retaining the term julq {d) instead of /xq 
in the first term of the right-hand side of the equation which may be written, on 
making use of (23), 



■X 



J (X) = Mo (0) Se-^°'H-x)sec?_ 1^^ J (X') .Ei { -Ko (X-X')} dX' 



+ J{X')Ei{-K,{X'-X)}dX' . . (32) 

Jx 
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In the notation of the preceding section we have to solve the above integral 

equation for 

M (X) = J (X)/Smo (0). 

We have 

while 

Ei{-Ko{X-X')}dX'+ Ei{-K,{X'-X)}dX' 



(p (X) = — |/Co 







If we notice that 



/•CO 



M{-ax)dx=f{ac)/a, . . . . 



where 



"CO 



f{x)^ = e~''-\-xEi{-'x) = x\ e~%r'^du, . . 



• t < 



• • (34) 



(35) 



X 



we find 



^(X) = i[2-/(KoX)-/{Ko(H-X)}],co/Ko. . 



. . (36) 



This expression is symmetrical with respect to the plane X = ^H, where it has its 
maximum value. The minimum value of the expression occurs at X = and X = H. 
We write for brevity 






y- 



(37) 



We then find 



B= {l^/(iC)}c/C and 6 = i{l-/(C)}c/C. .... (38) 



We also have 



a = H"^ r e-^o^^-^^-^'^dX = (l-e"^-^0/C sec ^, 

Jo 



(39) 



or, introducing the notation 



we have 



G{x) = {l-e--')/x, 

a = G (C sec 0. 



.... 



« e 



»••»•• 



• . 



(40) 
(41) 



Further, we find 



H 





^ = H-^ <p (X) dX 



or 



O — "o XI 



{l-f(KoX)}dX+\ [l-/{Ko(H-X)}]c?X c/C, 

Jo J 



/3 = C-i 



■c 



{l^/(^)}^^ c/C. 



* The function / (x), as defined above, occurs in a number of absorption problems. Its properties are 
described in a paper by the writer (King, L. V., * Phil. Mag./ February, 1912, p, 245), where a short 
riumerical table of the function is given* 

yOL, CCXII. — A. 3 D 
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It can easily be shown by integrating by parts that 



so that 



Jo 



/3={l-i/(C)-iG(C)}c/C. 



. . . . (42) 



(43) 



9 3 » » S S » 



Finally we obtain for eg, e, €i the expressions 



^ 



'2 



y/C+/(iC)c/C 



G (C 8ec 



r/C+i{/(C)+G(0)}c/C 



, y 



ft 9 



e 



-Osect 



€i = 



y/C+|- {1 + f{0)} c/G 



J 



The approximate solution of (32) may now be written 



. . . (44) 






= e 



~Ko(H-X)sec t 



Mo 



+ j e I ^ (X), . 



. . . (45) 



where <p (X) is defined by (36) and the yalnes e^, e^ eg are given by (44) and are 
employed in (45^^ according as we wish to make use of the extreme or mean solutions. 
From (20) the radiation scattered to a point on the earth's surface^ contained in a 
small solid angle w in a direction (/> with the vertical, is given by 



(lO 



T 



(JO 



l?o 







J (X) e"^^"^ dK 



Since R = X sec ^S this equation becomes 



T (^, = sec cj> 



rii 



e 



-KoX8eo<i^ 



J (X) dX, (46) 







where the dependence of the sky radiation per unit solid angle on the direction ^ and 
on the zenith distance of the sun ^ is denoted by T (^^ ^). 

If we denote by ft)R(^5 f) the radiation to a point at X = H contained in a small 
solid angle w pointing earthwards in a direction (p with the vertical, we have 

R {<p, = sec ^ r e-Ko(n--x)sec<^ j (x) dX. ..... . (47) 

Jo 

Substituting for J (X) from (45) we obtain from (46) 



T (^, = Mo (O) S sec ^ 



e 



— KoH sec t 



H 



vO 



g-KoX(sec^-sec.d ^X+ ( € 1 J 6"^^^^^^^ ^ (X) dX [> ' 



V 



%. 



w 
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The first term in the brackets may be written 

g-Cseo* Q |(geg ^_ge(, ^) c} C/Ko if ^ > </. 

or 

g-cseof Q {{sec ^_sec ^) 0} C/Ko if ^ > f . 

The evaluation of the second integral is more difficult ; we have 



**•«••• 



(48) 



J^g-KoXsec*^(X)(^X=igj e-«»^-<'*[2-/(KoX)-/{Ko(H-X)}]c^X. . (49) 



Writing KqX = u, the above integral takes the form 



1 c 



ZJ\.q \J 



-2(l-,g-Csec<f>) 

sec 9^ 



■c 



e' 



■u sec 



'^/(t^c^t^-e-^^^^* f e''''''^f{u)du . . (50) 

Jo 



We denote by B(aj) and B(— a;) the functions 

B {x) = -S^* (— 0?)— log aj, B {—x) == -Sa^ (cc)— log X, . 
The expansion when x is small of the exponential integral is 



■ (51) 



Ei{x) = y-^-flogx^-bx-h ^^-^i + 



I . , t I " • • • « » • 



2.2! 3.3! 4.4! 



. (52) 



where y is Euler's constant, y = 5572 and the expansion holds for both positive and 
negative values of x. 

Thus when x is small the expansions for B (x) and B{—x) are 



'B{x) = y—x-^ 



X' 



X' 



and 



2.2! 3.3! 






X' 



B(-~i:c) = y-^X-i r + 



X' 



2.2! 3.3! 






> 



J 



« • « • 



(53) 



It can be shown that the integral 



1 



+ 



a 



1 



log 



a 



a-\-h 



^e^^'Ei{-ac)+j^Ei{-{a-\-h)c) . . (54) 



The result holds for positive as well as negative values of h provided the argument 
of (a + 6) c in the logarithm and in Ei { — {a-\-h)c} is the same. 

By means of this result and in the notation defined in (51) it can be shown that 
(50) reduces to 

. . (55) 



x2- Ji_ __JL_<T) (n rhS 



« • e • 
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where 



(l_e-csec.){i_j(c)| 

+ cos^[B(C)-B{C(l+sec^)}+e-*^-'=*{B(C)-B[-C(sec^-l)]}].. (56) 



We notice that for ^ = 0, see <p = 1, 

$(C,0) = (l-e-°){l-/(G)}+B(C)-B(2C) + e-«{B(C)-y 
where y here stands for Euler's constant. 
For (f) — ^TT, sec ^ = oo (56) reduces to 

$(C,l7r) = l~/(C) 



. (57) 



(58) 



From (48) and (55) the expression for the intensity at the earth's surface of 
radiation from that portion of sky which is in any direction ^ with the vertical and 
azimuth i//> measured from a vertical plane through the sun is, for ^ >^. 



T(^,0 = 



Mo(^) 

JlVa 



s 



Csec^e-°^'=°*G{C(secf-sec9i)}+ig( e )$(C,^) 



,62. 



(59) 



while for f <^ we have 



T i'P, 



Mo 



(0) 



K, 



C sec <pe-'^''''<= G {G (sec ^-sec 0} +1 




$ (C, <!>) 



(60) 



From equations (3) and (7) we notice that 



so that the factor 



Mo (^) = I (1 + cos^ ^) Mo = i (1 + cos^ &) fo/^TT, 

SMo(0)/Ko = S(4,r)-^f(l+COS^0)c/C. . 



B • » 



, (61) 



From the polar diagram fig. 1 we see that 



cos = cos ^ cos f + sin <p sin ^QOS\fr, (62) 

The first terms in equations (59) and (60) give the contribution to the intensity of 
sky light due to the sun s radiation which has been once scattered by the atmosphere. 

If G is small and if the attenuation is due to scattering only, (c = C), the first 
term of (60) gives for the scattered radiation coming from a direction ^ the value 



(47r)-^S 1(1+ cos^ e) c sec 9^ . 



(63) 



which agrees with the value obtained by Kelvin/^ 

The first terms in (59) and (60) taken as they stand with C = c take into account 
the fact that both the incident and scattered radiations suffer attenuation : this effect 



■^ Kelvin, loc, cM, equation (19)^ p. 313. 
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is considered in a formula obtained by Raylbigh/^^ The term y in C = c + y allows 
for attenuation by absorption alone {i.e., without scattering). The second terms in 
(59) and (60) represent the contribution of self-illumination to the scattered radiation 
coming from any particular direction. An evaluation of this effect has not, so far as 
the writer is aware, been submitted to calculation although the importance of the 
effect is realized both by Kelvin t and Rayleigh,"^ and in an analogous problem by 

LOMMEL.J 

The expressions (59) and (6O) for the scattered radiation from any direction besides 
depending on the coefficients of absorption for the radiation of wave-length under 
consideration, depend also on the angular co-ordinates of direction ^ and \fr as well as 
on ^ the zenith distance of the sun. If we consider the intensity from zenith sky the 
expressions are greatly simplified. Writing </) = and = ^, (59) gives, 



T(0,0 = ^S 



c 




C6-^G{C(sec^-l)}+i--^( 6 )*(C,0) , . . (64) 



where $ (C, O) is given by (57) and is tabulated in Table V. 

It is not difficult to construct a double-entry table in terms of C and ^ giving the 
values for the functions which occur in (64) so that observations on zenith sky are 
most appropriate for comparison with the results of calculation from the attenuation 
coefficients determined by observations at the same time. 



TT 




The intensity of sky radiation from the direction of the horizon ( = - ) is given by 

T:(Tj'j = ^^ie-'^-^+^^(l']^{G,h)[, .... (65) 

where by (62) 6 is given by cos 6 = sin f cos ^ and $(C,-|7r) is defined by (58). It 
will be noticed that this formula unlike the approximate one (63), obtained by Kelvin, 

* Eayleigh, ^Phil. Mag.,' XLL, p. 116. 

t Kelvin, loc, cit,, p. 302, sect. 54. 

J LoMMEL, E., 'Sitzungb. der math.-phys. Glasse der K. Bayer. Akad. der Wiss.,' Bd. 17 (1887), 
p. 95; (the analysis is reproduced by Muller, ' Photometrie der Gestirne,' Leipzig, 1897, pp. 47-52). 
In working out the scattering of radiation by opaque, diffusely reflecting surfaces on an absorption 
theory, Lommel obtains a first order approximation to the efl'ect of seK-illumination in giving rise to 
deviations from Lambert's Law of diffuse reflection, the modified formula being referred to as the 
Lommel-Seeliger Law of Illumination. This problem is included as a particular case of the investigation 
of the present paper. Lommel's formula for diffuse reflection is represented in the present instance by 
equation (69), which is the solution of an integral equation expressing exactly the effect of self -illumination. 
MtJLLER points out that even Lommel's modification of LamIbert's Law does not represent exactly the 
results of observation on diffuse reflection ; it may happen that the more complete solution represented by 
(69), adapted to the case of intense absorption, gives a better representation of fact. This point must, 
however, be left over for further investigation. 
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remains finite, although its apphcation to the case of the earth's atmosphere is some- 
what invahdated by the curvature of the earth. 

On examining (59) it will be noticed that as long as G is not too large we may 
write 



T {i>, = MS sec ^T (0, = r^l sec ^T (O, 0, ■ 

MoUJ l+COS f 



(66) 



which may be taken as a rough approximation for the intensity of sky-radiation from 
any direction not too close to the horizon where (65) must be used. When C is large 
the wave-length corresponding is small, and the intensity in the normal solar spectrum 
outside the atmosphere also becomes small. Hence the formula (66) is sufficiently 
accurate when the total intensity on a horizontal plane is required. Denoting by 
H (f) the intensity of scattered radiation of wave-length X received per unit time on 
unit area of horizontal surface, we have 

H (f ) = T (^, f) cos0(iw, 



the integral being taken over the entire sky. 
Since dw = sin (f> d<l> dyjr^ we have 



H (f) = d^j/^ T (^, f) cos sin ^ dcj)^ . 

Jo Jo 

and making use of the approximation in (66) we have 



(67-) 



T(0,n 



'27r 



"VsT 



H (0 = ^ ^ ' Y, dyj^ (1 +cos' e) sin ct> dct>. 

l + cos f Jo Jo 

The total intensity for all wave-lengths per unit area of a horizontal surface is 
given by 

H(0c?X = 3, ^^ ,A T{OA)d\ . (68) 

f(l+cos^gJo 

Of some interest is the intensity of the radiation which is scattered from the 
atmosphere back into interplanetary space. If we write (H — X) for X in (47) we 
notice that 

E (0, ^) = sec <t> r e"KoXsecc|> j (H-^X) dH ; 

Jo 

we also notice from (36) that 0(H--X) = (X). 
We thus obtain from (45) 



R (0, t) = Mo (0) S sec 
which reduces to 



rll 

I g - Ko X (sec <|) + sec t) ^^ i 

Jo 




g--K„Xsee*^(X)c^X 



R (^, ^) = ^foMS 



K 







C sec ^G {C (sec ^ + sec 0) +i" A ( ^ ) * (C, ^) 



• (69) 
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We notice that R(0, ^ is only approximately equal to T {(p, f) when C is small, i,e.^ 
that the usual assumption that as much radiation is scattered in the direction of 
incidence as in the opposite direction only holds approximately when the coefficient of 
attenuation is small. 



§ 6. Note on the Polarization of Sky Radiation, 

It is well known that sky radiation is partly polarized in a vertical plane passing 
through the position of the sun (the principal plane) : in so far as the radiation to be 
scattered is direct solar radiation, the polarization ought to be complete. That 
portion of the sky radiation due to self-illumination is largely unpolarized and may to 
a large extent account for this deficiency from complete polarization : this point is 
mentioned by Rayleigh in his 1871 paper and the analysis of the present paper 
enables the magnitude of this factor to be roughly estimated. The complete solution 
of the problem from this aspect would require us to split up the incident radiation 
into two components, one of which is polarized in the principal plane, the other at 
right angles to it: the effect of self- illumination would lead to two simultaneous 
integral equations in three variables, the solution of which would be much too 
complicated to be useful. 

If, however, we refer to equations (59) and (60) it will be noticed that the 
expression for the intensity of sky radiation may be written in the form 

T{^,O = Mo(e)S{P(0,f) + Q(^,O}/Ko . ..... (70) 

where P {<p, C) stands for the first term in the brackets of (59) or (60) and 



2/ 



represents the effect of self-illumination. 

In default of a rigorous solution it is not unreasonable to suppose that the portion 
of the scattered radiation due to self-illumination is independent of the angle of 
polarization of the incident radiation. As far as the primary scattered radiation alone 
is concerned, the intensities polarized in the principal plane and in a plane at right 
angles to it are in the proportion 1 to cos^ 6. Thus from (70) the ratio of the sky 
intensities polarized in the principal plane and in the plane at right angles to it are 
given by the ratio 

If we make use of the approximations of equation (66), 

Tif,0 = sec <p P (0, Mo {0)/h iO, Q{i>, = sec f Q (O, () mo (0)/mo (() 
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(7l) may be written 

T. {<!>, _ cos''0 + Q(O,O/P(O,O ,^.^ 

T,(^,0~ 1 + Q(0,0/P(0,0 ^^ 

From this formula we see that sky polarization is most complete for 6 = |-7r, ix,, 
over a great circle an angular distance |-7r from the sun polarization would be 
complete in light from this direction if Q (O, were zero, i.e., if the effects of self- 
illumination were negligible. 

For the sake of comparison with observation the ratio T^ (^, ^/Tg {^, C) is calculated 
for zenith sky, {(p = O), in the form 

T, (0, _ i+Q(o,0/P(o,0 (^.. 

%{0,0 cos^^+Q(0,0/P(0,0 ■ ^'''^ 



Papt TTT 

§7. Analysis of Ohservatiofis on the Attenuation of Solar Radiation by the Earth's 

Atmosphere. 

The intensity of solar radiation transmitted to a station at a height x above sea- 
level is given by equation (23) in the form 

E(X) = S6-^«^^^""^^«^^^, 
where 






/'CO 



^ dx, and H ^ ^ dx 



po 



Po 



represents the height of the " homogeneous atmosphere/' 
Thus 



""".xx. ^^^ ax, so Tinau 



>G0 

pdx 

X 



p 



Jo 



where p and po ^^'^ the pressures of the atmosphere at the station and at the sea-level 
respectively. We thus have, writing C^, = Cp/po, the expression 

E(ir)-Se-^^^^«^, (74) 

If the heights of the barometer at the station and at sea-level are known at the 
time of observation, the comparison of absorption coefficients is independent of the 
law of variation of atmospheric pressure-gradient with height. 
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We notice from (37), writing C = c + y = *:„H + ao]iIo, that 



t, (y, J 



> . . . . (75) 



which gives 



^G - ^X-' + y , where /3 = ^^^ {n, - 1 f H/No, J 



Qt.Po/p = fiX^^ + y (76) 



Extensive observations on the determination of the coefficient of attenuation C^, for 
different wave-lengths have been carried out at various stations by the work of the 
Smithsonian Astrophysical Observatory."!' The mean coefficients of atmospheric 
transmission for Washington, Mount Wilson, and Mount Whitney have recently been 
given by Abbot.I These are quoted in Table I., while for the sake of independent 
comparison, the results of MtJLLER§ for Potsdam are also added. 

In order to study the correctness of formula (76) the coefficients of transmission for 
the different stations are plotted on a base X~^ X being measured in microns (lO~^ cm.). 
The results are shown in Diagram I., and give rise to a number of straight Unes. If the 
absorption were due to air-molecules alone, (76) shows that we should obtain a family 
of straight lines all passing through the same point (X"^ == -y/^). The straight lines 
actually obtained show that some variable factor in the atmosphere other than the 
molecules themselves is effective in attenuation, especially for stations below the level 
of Mount Wilson. This factor is generally referred to as atmospheric '' dust.'' A 
slight generalization of the analysis by which (76) was obtained enables us to interpret 
the results shown graphically in Diagram I. 

Let N^ be the number of " dusf particles per unit volume at a height x above the 
earth's surface. The coefficient of attenuation, which includes the effects of scattering 
and absorption both by air molecules and by '' dust," may be written in the form 

C\, -(/9X~'^ + y)|)/po + K^o(x) ^dx, (77) 

where N^o is a constant representing the number of dust particles per unit volume at 
the earth's surface, and K^o M depends on the nature of absorption and scattering by 
the dust particles. The distribution of '' dust " N^ in the atmosphere may be regarded 

* The representation of the coefficient of attenuation as the sum of two terms, one constant and the 
other varying inversely as the fourth power of the wave-length, seems to have been first recognized by 
Becker (see Kelvin, 'Baltimore Lectures,' p. 321, equation 33) from an analysis of Mijller's 
observations. 

t ' Annals of the Smithsonian Astrophysical Observatory,' vol. IL, by C. G. Abbot and F. E. Fowi.E, 
Washington, 1908. Eef erred to subsequently as ' Annals/ vol. IL 

I Abbot, C. G., "The Sun's Energy-Spectrum and Temperature," ' AstrophysicalJournal,' XXXIV., 
October, 1911. 

§ MtJLLER, G., " Photometric der Gestirne," Leipzig, 1897, p. 140. 
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as an unknown and variable quantity. Diagram I. shows th at for Mount Wilson and for 
Mount Whitney the effect of " dust '' is small, i.e,^ we may take N' = for x greater 
than 1780 metres. In this case i^7) reduces to (76) and the attenuation may be 
taken to be due almost entirely to the effect of air-molecules, while the existence 
of a small term, y, indicates that eyen in the comparatively dust-free air above Mount 
Whitney there is a small amount of attenuation by absorption, i.e.^ a direct conversion 
of solar radiation into thermal agitation of atmospheric molecules. 

We notice from Table II. that for the comparatively dust-free air above the 
level of Mount Wilson the value of y under standard conditions of pressure is 
y = ao H = *032, or, since H == 7*988 x 10^ cm. at 0" C, the value of a^ is 

ao = 4-Ox 10-""^ cm."-\ (78) 

On referring to (8) we notice that a^ is the fraction of radiant energy converted per 
centimetre of path into thermal molecular agitation. This fraction is greatly mcreased 
by the presence of small solid particles such as '' dust," &c. From (9) we can estimate 
the rate of increase of temperature in a gas under atmospheric pressure due to solar 
radiation passing through it. 

We have -~:^ = "^- from equation (9). 

at pqS 

Taking the value of a^ for dust-free air from (78) and writing E = 1*92 calories per 
minute, p^ = '001293, 5 = '237, we find 



-J- = 2*5 X 10""^ degrees (1 per minute 
= 1*5 X 10"^ degrees CX per hour 



>, 



, (79) 



In the case of ordinary air at sea-level, the Washington observations from Table II. 
show that the rate of increase of temperature calculated in (79) must be increased to 
about six times this value. 

With regard to sea-level stations, Diagram I. seems to indicate, both for the 
Washington and Potsdam observations, a marked change in the nature of the 
absorption due to ''dust'' in the neighbourhood of 'GlO/x. We therefore discuss 
separately the case of long- wave radiation (X>*filO/x) and short-wave radiation 
(X<*610m). 

(i) Long-wave Radiation (X^'GlOfx), 

For long waves the straight lines of Diagram I. show that the term Wq (x) of (77) 
must be of the form 

where /3^^ and y^ are constants for the range of wave-lengths greater than *610/x. This 
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indicates that the presence of ^' dust" gives rise to both absorption and scattering ; 
{^7) may then be written 



where 






1 — y ~Y- — y 



(ajQC* 



JP ^ X iN 



Diagram 1. shows clearly that for long waves the straight lines both for Potsdam 
and Washington intersect in the same point, X"'* = — y//8, as do the lines for Mount 
Wilson and Mount Whitney. 

From (80) this result requires that 

jjjy z=z ^ jy ^ , . . , . . . . . , (81) 

a condition which is independent of the law of distribution of the " dusf particles. 
The ratio ^ is proportional to the ratio 

r 

energy of incident wave scattered by small particles 
energy of incident wave converted into molecular agitation ' 

and (81) indicates from the results of observation that for long-wave radiation this 
ratio is independent of the nature of the scattering particles, whether '^ dust " or air- 
molecules. This result throws some light on the question raised at the end of §1 as 
to the mechanism by which the molecules of a gas can convert a portion of the 
radiant energy incident upon t-hem into thermal molecular agitation. The same 
mechanism which is effective in scattering radiation is also capable of effecting 
molecular velocities, and therefore the rate of increase of temperature in such a way 
that the ratio /3/y is independent of the nature of the molecule or even of the small 
" dust ^^ particle giving rise to the absorption and scattering. 

(ii) Short-ivave Radiation (X<"610m). 

For short waves we may suppose that the incident radiation is not scattered by 
the dust-particles but is absorbed and converted into heat. On this supposition K^^ (^) 
in {77) is of the form y^, where y is a constant independent of the wave-length. 
Equation {^7) then takes the form 

C, = X-^/3p/po + rPM) + r :^dx = X-"'^p/p,-^rp/p,. . . . (82) 

The slope of this line in the graphical representation of Diagram I. is given by 

tan 0^ = ^p/Pq .... ...... (83) 

and is independent of the distribution of " dust " in the atmosphere. This conclusion 
is justified by calculating /3 for the various stations from the slope of the lines in 

3 E 2 
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Diagram 1. and a knowledge of the mean barometric pressures at these stations. The 
results are tabulated and described in Table IT From these values of ^ we may by 
(75) calculate N^ the number of molecules in a gas under standard conditions of 
pressure and temperature. We take r\^—l ^^ °00()293, H ^ 7 '988 x 10'' cm. at 0° C, 
and (71) then gives 

N, -^ 2*269 xlO^V/^- ' • (8^1) 

The results for each set of observations are given in Table 11. The values of N^ agree 
remarkably well among themselves, and give for a mean value'^'^ 

This result is in substantial agreement with the value of Nq obtained by 

Rutherford and GEiOERt 

No -= 272 X 10^-', 

and with the value obtained by MiLLrKAN| from a recent determination of the 
elementary electrical charge, e = 4*891 x 10"^ E.S.U,, which gives 

No = 2*644 X 10'', 

taking the Faraday constant to be 9*655 absolute E.M. units. Tliis agreement 
indicates that for short-wave radiation the scattering is almost entirely due to air- 
molecules, while the effect of "dust'' is to produce a genuine absorption effect, i,e. a 
direct conversion of radiant energy into heat. 
The term 

:z3 y -f- ^—^ y ——- (IX . . , . . . . . . (86) 



p J. 



i) 



involves the distribution of atmospheric "^ dust' and may therefore be expected to be 
an extremely variable factor with respect both to place and time. The method of 
analysis of the present section offers a convenient method of studying the variations 
in the distribution of atmospheric " dust " and their connection with other 
meteorological phenomena. 

* The close agreement of coefficients of attenuation calcalated from tlie formula c = Pj^"'^^ using 
EuTHERFORD and Geiger's value of No, with coefficients calculated from observations on selected clear days 
at Washington and Mount Wilson, was pointed out by Schuster (" Molecular Scattering and Atmos- 
pheric Absorption," 'Nature,' July 22, 1909; 'Optics/ 2nd ed., 1909, p. 329). 

t EuTHEREORD, E., and Geiger, II, " Gliarge and Nature of the a-Particle," ' Koy, Soc. Proc.,' A, 
vol 81, 1908, p. 171. 

I MiLLlKAN, E, A., " The Isolation of an Ion, a Precision Measurement of its Charge, and the 
Correction of Stokes' Law," ' Phys. Eev.,' XXXII., April, 1911. A summary of the various physical 
measurements which lead to the value of No is given by Perrin, 3. (' Annales de Chimie et de Physique,' 
8"^^ ser,, September, 1909, translated by Soddy, F. ; ' Brownian Movenient and Molecular Reality,' Taylor 
and Francis, 1910, p. 90). 
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§ 8. On the Intensity of Sky Radiation as Calculated from the Mean Coefficients of 

Attenuation at Mount Wilson and Washington, 

It will be noticed from (64) that the intensity from zenith sky for different wave- 
lengths can be expressed in terms of the zenith distance of the sun and the coefficients 
of attenuation C, c, and y which are determined for a given station by an analysis of 
atmospheric transmission observations according to the method of the preceding 
section. In order to simpHfy the calculations from (64) the various functions of 
C and f which occur in this formula are tabulated in Tables III., IV., and V. 

The intensity from zenith sky is then worked out for the two stations Mount 
Wilson and Washington from, the mean coefficients of attenuation at these two places. 
In the first case the numerical values corresponding to the extreme and mean solutions 
of the integral equation are carried throughout all the calculations and in this way 
give the limit of errors due to an approximate solution. In the second case, which is 
taken as typical of a sea-level station, numerical values corresponding to the mean 
solution of the integral equation are alone given, the reason being that the absorption 
at sea-level is an extremely variable quantity which would give rise to fluctuations 
in sky radiation probably exceeding the difference of the extreme solutions. 

Zenith intensities of sky radiation for different wave-lengths and for vsirious zenith 
distances of the sun are given in Table VI. for Mount Wilson and in Table VIII. for 
Washington. The unit of intensities is arbitrary and is that employed in Table I. for 
the normal solar « spectrum outside the earth's atmosphere. The results are shown 
graphically in Diagrams II. to VII. for Mount Wilson and in Diagram XII. 
for Washington. From these curves it is possible to obtain by double interpolation 
the intensity from zenith sky corresponding to any wave-length and zenith 
distance of the sun. In this way a comparison was made of the quality of sky 
radiation obtained by calculation v/ith that obtained experimentally at Mount 
Wilson.'^ The results are given numerically in Table X. and are compared 
graphically in Diagram XIII. 

By integrating the curves for the zenith intensity and making use of the 
approximate formula (66) we are able to obtain a rough approximation to the total 
intensity of sky radiation from any direction. On comparing (66) with (65), for which 
(j) = ^TTj it must be noted that the first of these formulae can only be used for zenith 
distances which are not too great. It will also be noticed that T {^tt, does not 
vary as rapidly with the wave-length as T (O, (*),. The nature of the scattered 
radiation from a portion of the sky near the horizon is governed principally 
by the term Se"^'''''^, i,e,, it is of nearly the same quality as the direct sunlight 
except for the contribution to the intensity due to the second term in T (|-7r, ^) 
which represents the effect of self-illumination. This feature is not shown by 



■K- i 



Annals,' vol. IL, Table 32, p. 155. 
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expressions for sky radiation hitherto obtained ; in fact tlie whitish colour of the 
sky near the horizon was considered by KEnviN"^ to be an objection to the theory of 
scattering. 

By calculating the values of the intensities of solar radiation which reaches the 
earth's surface for various zenith distances of the sun and for various wave-lengths, 
and by integrating the curves so obtained, we obtain the values of the total intensity 
of solar radiation reaching the two stations Mount Wilson and Washington, Taking 
the value of the solar constant to be 1*922 calories per square centimetre per 
minute, t the results are given in the same units and are calculated for solar 
radiation incident on a plane normal to the sun's rays as well as for the solar radiation 
incident on a horizontal plane. The results are given in Tables VII. and IX., and 
are shown graphically in Diagrams IX. and X. 

These results enable the value of the total sky intensity in any direction to be 
compared with the intensity of direct solar radiation at the station in the form 
sky /sun for equal solid angles (semi-diameter of sun taken as !(/ of arc). The results 
are given in Tables VII, and IX., both for Mount Wilson and Washington, and are 
shown graphically in Diagram XI, 

By interpolation from the diagram just mentioned the values of sky/sun for 
various directions of sky and for various zenith distances of the sun are compared 
with the values observed at Mount Wilson.| The comparison is made in Tables XI. 
and XIII. ; the results of calculation are in fair agreement with observation 
except for regions of sky near the sun and near the horizon. The first of these 
discrepancies is probably brought about by the simplifying assumption made in 

equation (l7) by writing jl for /x (0) and iuL{rr^). The existence of such a term 
depending on an angular co-ordinate complicates the integral equation beyond hope 
of solution ; it can be seen, however, that its existence gives rise to a bright region 
of sky in the neighbourhood of the sun not represented by the approximate solution 
considered in the present paper. The discrepanc}/- which exists in the case of direc- 
tions of sky near the horizon is due to the failure of the approximate formula (66) 
for large values of (p. In such cases the complete formxula (60) should be employed. 
The use of this formula would require the tabulation of functions of two and three 
variables, which might be undertaken when more numerous and more accin^ate 
observations on sky radiation are available. It must be remembered in making these 
comparisons that the result of sky observations on certain specified days are com- 
pared with values based on mean coefiicients of attenuation. The only satisfactory 
method of making the comparison is to obtain sky radiation observations on the 
same day as observations are made for the coefficients of attenuation of solar 

■^ Kelvin, he. dt,, p. 307. 

t Abbot, 0. G., and Fowle, F. E., ''The Value of the Solar Constant of Eadiation," ' Astrophysical 
Journal,' XXXIIL, April, 1911, p. 191. 
+ 'Annals,' vol. IL, Table 32, p. 15L 
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radiation/^ In the present paper no account is taken of the reflectmg power of the 
earth's , surface. The tolerable agreement between the results of observation and 
those of a theory based on a non-reflecting surface shows that for an ordinary land- 
scape the effect of reflection from field and foliage need not be as great as is sometimes 
supposed. t The effect of snow on the polarization of sky radiation is well-knownj : 
on this point Lord Rayleigh mentions in his 1871 paper the interest which woifld 
be attached to sky radiation observations taken over a landscape' covered with snow 
or over the sea, reflection from these surfaces being in both cases especially 
determinate. § 

For meteorological purposes it is important to know the total solar radiation 
incident on a horizontal plane as well as the contribution due to sky radiation. The 
results are given in calories per square centimetre per minute for Mount Wilson and 
Washington in Tables VII. and IX., and are shown graphically in Diagrams IX. and X. 
It will be noticed that for large zenith distances of the sun the contribution of sky 
radiation to the total radiation on 'a horizontal plane is a very considerable fraction of 
that due to direct solar radiation. The factor of sky radiation is thus of considerable 
importance in the meteorology of northern latitudes. 

The agreement between the results of calculations and such observations as are 
available gives rise to a hope that the present communication may serve as a guide 
towards systematic observations of the type dealt with, and to their interpretation in 
terms of a theory of scattering and absorption ; by this means one may hope to obtain 
absorption constants and methods of using them which will be of some service to 
meteorology and astrophysics. 

Summary, 



The analysis of the present paper seems to support the view that at levels above 
Mount Wilson molecular scattering is sufiicient to account completely both for 

■^ In the " Eeport on the Astrophysicai Observatory " (G. G. Abbot, * Annual Eeport of the Smithsonian 
Institution,' 1911, p. 65), the Director announces that sky radiation observations have been successfully 
taken at Mount Whitney (August, 1910). Since the transmission coefficients at the time of observation 
are also determined, the results will enable an accurate comparison of sky radiation results to be made 
with the values obtained by calculation in terms of the coefficients of attenuation. 

t NICOLS, "Theories of the Colour of the Sky," ^Phys. Eev.,' XXYL, June, 1908, p. 507. 

J Observations have been made by McConnel, J. C, on the effect, of the nature of the ground on the 
degree of polarization of the sky at 90 degrees from the sun ; the effect of a covering of snow is to 
diminish the degree of polarization. (Perntner, 'Meteorological Optics,' Part lY., 1910, p. G43. The 
above reference is taken from 'A History of the Cavendish Laboratory,' Longmans, Green & Co., 
1910, p. 129.) 

§ The analysis of the present paper could be extended without difficulty to take into account the effect 
of reflection from the earth's surface on the intensity and polarization of sky radiation, provided the 
landscape were covered by a layer of uniform substance such as snow, or observations were taken over the 
sea, which presents a determinate reflecting surface. 
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attenuation of solar radiation and for the intensity and quality of sky radiation. 
Even at sea-level the effect of '^ atmospheric dust " can be taken into account in a 
simple manner in the formulae for absorption and scattering. Should future observa- 
tions support the validity of the simple law expressed by equation (2) connecting the 
coefficient of attenuation with the wave-length, we may with considerable assurance 
make use of the law to obtain the coefficients for very short or very long wave-lengths 
when the direct method of calculation from high and low sun observations leaves room, 
for considerable uncertainty owing to the small intensities in the solar spectrum at 
these wave-lengths and owing to other experimental difficulties. 

Schuster'^ points out in this connection the extreme importance of determining 
accurately the form of the solar intensity curve outside the earth's atmosphere for 
short wave-lengths, since the effect of a solar atmosphere in absorbing and scattering 
radiation is to give rise to an intensity-curve which does not agree with that given by 
Planck's formula especially for short wave-lengths. 

Absorption and scattering of radiation by the sun's atmosphere, taken in conjunc- 
tion with effects of self-illumination, constitutes a problem analogous to that just 
considered for the earth's atmosphere. By making a comparison between the 
calculated variation of intensity of radiation of different wave-lengths over the solar 
disc and the results of observation it will be possible to determine from the intensity 
curve of the normal solar spectrum outside the earth's atmosphere the intensity-curve 
at the radiating layer of the sun. This corrected curve may then be compared with 
that given by Planck's formula and a closer approximation made to the temperature 
of the sun than the values now given. This investigation the writer hopes to be able 
to deal with in a future communication 



* Schuster, loc. cii 
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Numerical Tables and Comparison of Theory with Observation. 

Discussion of Results. 

Table I. — Coefficients of attenuation for Washington, Mount Wilson, Mount Whitney, and Potsdam. 
II. — Constants of atmospheric scattering and absorption. 
III.— -Tabulation of the function Ce~^ G {C (sec f -~ 1)}. 
IV. — -Tabulation of the function G (C sec (). 
V.^ — ^Tabulation of the auxiliary functions <E* (C,0), G(C), /(C). 
YI.— Relatiye intensities from zenith sky calculated from Moimt Wilson observations. 
VII.-— Total solar and sky radiation calculated from Mount Wilson observations. 
,, VIII. — Relative intensities from zenith sky calculated from Washington observations. 
IX. — Total solar and sky radiation calculated from Washington observations. 
X.— Quality of sky radiation at Mount Wilson, October 17, 1906. 
XL— Total sky radiation at Mount Wilson, October 19, 1906. 
„ XII.— Average intensities of sky radiation, Mount Wilson. 
„ XIII. —Polarization of sky radiation calculated from Washington and Mount Wilson observations. 
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Tablis I. — Coefficients of Attenuation for Washington^ 'Mount IVilson, Mount Whitney^ and Potsdam. 

The table of constants of solar intensities and coefficients of attenuation for Washington, Mount Wilson, 
and Mount Whitney is taken from the revised reductions based on the most recent observations of the 
Smithsonian Astrophysical Observatory."^ The unit expressing the relative intensities of solar radiation 
between wave-lengths A and X-^dX is arbitrary, and can be reduced to calories per square centimetre per 
minute by multiplying by a factor tj such that 

foe 
S c?A = 1 • 922 calories per square centimetre per minutest 

In terms of the unit of intensity given in the table, and taking dX = • Ifo, we find by graphical integration 

D 

S^a - 3321-0, so that ^ = -000578. 

No account is taken of the water-vapour bands which occur mainly in the infra-red, where the intensity 
of solar radiation is comparatively small. The position and extent of these bands is shown in the curve 
of the intensities in the normal solar spectrum given in Vol. 11. of the * Annals.*| 

The transmission coefficients for Potsdam are those given by Muller.§ 



* Abbot, 0. (jc., " The Sim's Energy-Specfcram and Tempera tore," 'Astrophysical Journal,' X.XXIY., October, 1911, p. 197. 
t Abbot, 0. 0., and Fowle, F. E., " The Value of the Solar Constant of Eadiation," ^Astrophysical Journal,' XXXIIL, 
April, 1911, p. 191. 

J * Annals,' vol. IT., p. 104. 

§ MtJLLEB, G-., 'Die Photometrie der Grestirne,' Leipzig, 1897, p. 138. 
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Table I 



Wave- 
length in 
microns, 
X. 



Mean 

solar 

intensity, 

S. 






'30 fi \ 





•325 ■ 





■35 





•375 





•39 





•42 





•43 





•45 





•47 





•50 





•55 





60 


0' 


70 





80 


1 


00 


1- 


30 


1- 


60 


2- 


00 


2- 


50 


3- 


00 



(440) 

1285 

2700 

3464 

3620 

5261 

5340 

6047 

6253 

6064 

5627 

5047 

3650 

2672 

1664 

897 

526 

245 

43 

(12) 



Per- 
centage 
probable 

error. 



30 



7 


•3 




, 


2 


•7 




1 


•7 





•445 


1 


•3 





•586 


1 


•7 





•600 


1 


•4 





•640 


1 


•8 





•671 


1 


•9 





•705 


2 


•1 





•739 


2 







'760 





•4 





■839 


1 


•2 





■865 





•7 





■901 





7 


0" 


916 


1' 


4 


0- 


930 


2- 


4 


0" 


909 


4- 


8 


0^ 


870 


45- 






— 



Mean coefficients of atmospheric transmission, e-^\ 



Washington, 
1902-07. 



e 



-Cr 



^iB' 



810 

534 

511 

446 

399 

350 

302 

274 

176 

145 

104 

0876 

0726 

0943 

139 



Mount Wilson, 
1909-10. 



e 


























(550) 

612 

662 

694 

764 

778 

800 

827 

858 

876 

890 

942 

964 

973 

972 

975 

957 



0^(900) 



G 



W' 



Mount Whitney, 
1909-10. 



g-cx 



Cn.. 



598) 
491 

365 

269 

251 

223 

190 

153 

132 

117 

0597 

0367 

0274 

0284 

0253 

0439 

105) 



(0 








































•510) ( 


•691) 


•584 


•538 


•660 


•416 


•738 


•304 


•763 


•270 


•806 


•216 


•822 


•196 


•851 


■161 


•880 


•128 


•900 


105 


•918 


0866 


•934 


0683 


•956 


0450 


•972 


0284 


•980 


0202 


•980 


0202 


•978 


0222 


•940 


0619 


•930 


0726 


•910 


0943 





Mean transmission coefficients for Potsdam. 




A. 


e 


• 


Go;. 


X. 


g-Ca-^ 


Git-. 
•186 




■44/x 





•706 


•348 


•58/y. 


0^830 




46 





■740 


•301 


•60 


0-840 


•174 




■48 





•764 


*269 


•62 . 


0-850 


•162 




•50 





■781 


•247 


•64 


0-861 


•150 




52 





795 


•229 


•66 


0-871 


■138 




54 


0- 


808 


•213 


•68 


0-881 


•127 


•56 ' 


0^819 


•200 









Coefficients of Jtienmition for Washington, Mount Wilson, Mount Whitney, and Potsdam. 
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Table IL— -Constants of Atmospheric Scattering and Absorption. 

The coefficients of attenuation given in Table I. for various stations were plotted on a large scale diagram 
against X"'^ as abscissae. The results are shown in Diagram I. From the large-scale drawing the slopes of 
the straight lines passing through the mean position of the observed points were obtained and the values of 
tan 6' and tan 6" corresponding to short- and long-wave radiation were calculated.* From a knowledge 
of the mean barometric pressures at these stations the values of p as defined in (75) were calculated for 
short-wave radiation, leading by (84) to an estimate of No, the number of molecules per cubic centimetre 
of a gas at 0° 0. and 760 mm. pressure. The values of the absorption constants F' and F" of (80) and (82), 
which include the effect of " dust " especially noticeable at low-level stations, are also given in the table. 
At levels higher than Mount Wilson the atmosphere is comparatively dust-free and F' ^ F" = y. This 
value of y can then be employed in (78) to give a numerical estimate of molecular absorption. 



[* If greater accuracy is required the lines of closest fit to the system of observed points can be drawn by calculating in 
each case the position of the major axis of inertia of the corresponding system of material points of equal weight. Formula} 
for the determination of this line are given by Kael Peaesof (' Phil. Mag./ vol. II., 6th series, E'ovember, 1901, p. 559), and 
also by Snow, E. C. ('Phil. Mag.,' March, 1911).— iVo/e added December 31, 1912.] 



Table HI.— Tabulation of the FmicMon^ Ce"*^G{0(sec{'- 1)}. 

The function G (x) = (1 - e~'')jx is tabulated by W. Lash Miller and T. K Eosebrugh^ in a set of 
extensive tables of the integrals 



'00 /• 3J roa 

,~-2 ^~-U r7n, n,-~l ^~U 



%i'^ e ^^ du, u ^ e ^^ du = ^i ( - x\ 



ue '^^d% u^e' '^Uin. 



Values are given to 9 significant figures at intervals of -001 between x ™ and x ^ 1, and at intervals 
of '01 between x = 1 and « = 2. From these tables the double-entry Table III. was easily constructed. 
The calculations of the present table, as well as those of tlie other tables, were performed on a slide-rule, 
so that their limit of accuracy is about one or Uvo parts in a thousand. 



* \¥, Lash Milldb and T. R. RosEBEuaH, ' Trans. Roy. Soc. of Canada,' 2nd serieB, vol. IX. , 1903, sect, iii., pp. 73-107 
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Tabt;e IL 



Station. 


Height 
above 
sea- 
level. 


Mean 
baro- 
metric 
pressure, 

mm. 

*752 
t763-7 
+ 617 
§446-7 


She 
tan^'. 


Potsdam .... 
Washington . . 
Mount Wilson . 
Mount Whitney 


metres. 

100 

10 

1780 

4420 


•00893 
•01016 
•00806 
•00592 



^. 



00904 
01011 
00994 
01003 



ition, X<-610jw. 


No. 


-105 
-192 
•032 
-032 


2-51x1019 
2-24x1019 
2-28x1019 
2-26x1019 



Long-wave radiation, 
X>-610p. 



tan (9". 


r". 


f% 


- r /r". 


-0170 


-054 




-32 


-0277 


-088 




-32 


-00806 


-032 




•32 


•00592 


-032 




•32 



Constants of Atmospheric Scattering and Absorption, 



Authorities :- 



* The mean barometric pressure at Potsdam is given by Mullee, Ioc. cif., p. 188. 

t The mean annual barometric pressure at Washington is taken from Barthoiomevv's 'Atlas of Meteorology/ 1899. 

X The mean barometric pressure at Mount Wilson is not given explicitl}^ in the 'Annals.' The value given in the above 
table is obtained by finding the reduction of pressure to sea-level corresponding to an elevation of 5886 feet (the air- 
temperature being taken at 60° F.) from Hazen's Tables. (' Professional Papers of the Signal Service,' No. YI., 
Washington, 1882.) 

§ The barometric pressure at Mount Whitney is obtained from observations given by Lang^ley, September 2-6, 1881 
(Langley, '^ Researches on Solar Heat : A Keport of the Mount Whitney Expedition," ' Professional Papers of tlie 
Signal Service,' 'No. XY.,1884.) 



Table III. 



c.\ 


! 

o^ 1 

1 


20". 


40°. 


60°. 


70\ 


80°. 


0-00 


•0000 


•0000 


• 0000 


0000 


■ 0000 


•0000 




05 


0476 


•0475 1 


0472 


0464 


•0443 


•0424 




•10 


0905 


•0902 i 


•0891 


•0861 


•0813 


•0720 




•14 


1217 


•1212 


•1191 


•1133 


•1072 


•0889 




18 


1502 


•1492 i 


•1461 


1377 


•1270 


1009 




22 


1766 


•1753 j 


1708 


1583 


1439 


1093 




26 


2002 


1985 


•1927 


•1761 


•1575 


•1150 




30 


2220 


•2200 


•2136 


1953 


•1749 


•1174 




34 


2420 


2395 


2298 


2052 


1776 


1198 




38 


2600 1 


2568 


2450 


•2163 \ 


1844 


1206 




50 


303 1 


•298 


•280 


238 ' 


1964 


1154 




60 


329 


323 i 


301 


248 \ 


1950 


1090 




70 


348 


340 


313 


250 


1930 


1006 




80 


359 


350 


318 


247 


1850 


0923 


-90 


366 


356 

1 


320 


241 


1740 


0856 



Table of the Function Qe ^ G {C (sec f - 1)} 
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Table 1Y .-^Tabulation of the Function G (C sec {"). 

This function occurs m calculating the mean coefficient, e, in the solution of the integral equation (44). 
It is easily tabulated from Lash Miller's tables. 



Table Y .—Tabulation of the Auxiliary Functions ^ (C, 0), G (G)^ /(G). 
The function /(.T.) defined by the relation 

u~^e~"^ du — e~^ + xEi ( - x) 

occurs in the theory of absorption in a flat plate which itself absorbs radiation from a uniform distri- 
bution of radiating elements throughout its own volume. A short table of the function f{x) has been 
given by the writer in a previous paper. ■^* 
We write as in (51), 

B {x) ^ Ei{- x) - log x^ B ( - a;) = Ei {x) - log x. 

The function B {x) is given in Lash Miller's Tables. f The function B ( - x) was calculated by making 
use of the tables of the exponential integral Ei {x). J 
From these a table of the function 

#(G,0) = (l-6^-^){l-/(C)}+B(C)-B(2C) + 6-^{B(C)-7} 
was constructed. 



* KiNa, L. V.5 " Absorption Problems in Eadioactivifcy," ' PbiL Mag./ February, 1912, p. 245. 
t W. Lash Milleb and T. R. Eosebrtj0h, loe.'cit., p. 81. 






Cf. Dale, ' Five-Figure Tables of Mathematical Functions/ Arnolds, London, 1908. 
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Table IV. 



X 

C. -x^ 


0\ 


20°. 


40\ 


60°. 


70^ 


80^ 


0-00 1 


' 0000 1 


• 0000 1 


• 0000 1 


• 0000 1 ' 


•0000 1" 


0000 




•05 


■975 


•974 


•968 


952 


'929 


873 




'10 


•952 


•948 


'940 


• 906 


•868 


770 




'14 


•933 


•929 


•913 


872 


'820 


688 




•18 


•916 


•910 


•891 


840 


•778 


•623 




•22 


•897 


•892 


•869 


809 


737 


•567 




•26 


•880 


•874 


•848 


•780 


•700 


•536 




30 


•864 


•856 


•828 


752- 


•666 


•476 




34 


848 


•839 


•808 


726 


633 


438 




38 


832 


•822 


788 


701 


•603 


•406 




50 


788 


777 


•735 


632 


•526 


•328 




60 


751 


740 


694 


583 


•472 


•281 




70 


718 


704 


656 


538 


•428 


•243 




80 


688 


672 


622 


499 


•386 


•215 


•90 


658 

1 


642 


588 


464 


•357 


•192 



Table of the Function G (C sec f). 



Table V. 



c. * 


(C, 0). 


G (G). J 


^(0). |{G(C 


')+/(C)}. HI 


+/(C)}. / 


©■ 


0^00 


•0000 


1 • 0000 1 


■0000 1 


•000 1 


•000 1 


•000 




•05 


•0095 


•9754 


•828 


•902 


•9U 


•897 




•10 


•0311 


•9516 


•722 


•837 


•861 


•828 




•14 


•0536 


•933 


•656 


•794 


•828 


•782 




•18 


•0801 


•916 


•600 


•758 


•800 


•741 




•22 


•1081 , 


•897 


•550 


•724 


•775 


705 




•26 


•1403' 


•880 


508 


•694 


•754 


672 




•30 


•1723 


•864 


469 


•667 


735 


641 




34 


2065 


•848 


435 


•642 


718 


613 




•38 


2414 


•832 


404 


•618 


702 


587 




50 


349 


•787 


327 


557 


664 


518 




60 


438 


•752 


276 


514 


638 


469 




70 


525 


•719 


235 


477 


618 


427 




80 


612 


•688 


201 


445 


601 


389 


'90 


686 


•659 


172 


416 


586 


356 



Table of Auxiliary Functions * (C, 0), G (C), /(C). 
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Table Y1,~— Relative Intensities from Zenith Shy Calculated from Mount Wilson Observations in Terms of Mean 

Attemiation Coeffidenis. 

The wave-lengths at Mount Wilson level corresponding to the tabulated values of the attenuation 
coefficients C were determined from the formula (76), = /3k"'^ + y, making use of the constants ^ and y 
given in Table 11. , which themselves are derived from observations on the transmission of solar-radiation. 

A large scale chart of the intensities in the normal solar spectrum outside the atmosphere was prepared 
from the data reproduced in Table I. The values of S corresponding to various values of A were then 
estimated from this curve of solar intensities. Auxiliary tables of the coefficients €25 e, e^ were calculated 
from (44), making use of Table V. In this way a table of the function 






ijM*^^^^) 



was prepared. This table, with Table III., enabled the term in square brackets in the expression for 
T (0, {■) to be calculated. Finally, making use of the values of S just determined, the present table of the 
intensity from zenith sky, T (0, {"), was calculated. The extreme and mean solutions of the integral 
equation are retained throughout. It will be noticed that the extreme solutions diverge rapidly for large 
values of and ( ; it must be remembered that for large values of C the value of S decreases with extreme 
rapidity, no appreciable intensity having been measured for a wave-length less than -S/x. Consequently, in 
calculating the total intensity of sky radiation on a horizontal surface, the divergence of the extreme 
solutions does not lead to a very great divergence in the final results. The results given in the present 
table are shown graphically in Diagrams II.-VII. 
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Table VI. 



c. 


X. 


S. 


0°. 


20°. 


40°. 


60°. 


70°. 


80°. 


•05 


1 

: 

•741/x 


3140 


9-8 
9^8 
9^8 


9-2 
9^2 
92 


7-7 
1-1 
7-7 


6-0 
6-0 
6-0 


5-1 
b-l 
51 


4-6 
46 


•10 


•571 
•510 


5380 


SO'9 
bO'6 
bO'2 


Jf.7-7 
47^4 
47^0 


89-8 
39^5 
390 


30-4 
30-1 
29-6, 


25-9 
25-4: 
248 


21-6 
20-8 
20-0 


•14 


6000 


88-5 
S7'2 
85^9 


83-0 
81-7 
803 


69-0 
^1-7 
664 


52-2 
50-7 
49-4 


U'6 
42^8 
413 


35-6 
32^9 
312 


•18 


•474 


6260 


118 •a 


116-7 
113^8 
1102 


96-5 
93-1 
90-7 


72-8 
69-8 
66-1 


61-2 
b1-b 
540 


48-0 

42^3 
381 


•22 


•446 


5910 


H6'0 
141^5 
136 8 


137-0 
132-5 
127 8 


113-0 
108-4 
104 


84-5 
79-8 
75-1 


70-1 
64^5 
59-6 


53-9 
45-4 
400 


•26 


•427 


5320 


166^9 
159-0 
151 


156-0 
148^5 
140 8 


128-5 
121-2 
114 1 


95-5 
81-1 
80-6 


79-4 
10-1 

625 


60-7 
4:8-0 
394 


•30 


•410 


4480 


168-0 
158-4 
148-5 


157-1 
147^9 
137^9 


130-0 
120^6 
108 2 


97-0 
86^6 
780 


80-5 
69^0 
603 


60-0 

4:3-1 

34-6 


•34 


•397 


4000 


173-2 
161^0 
148^1 


162-0 
150^0 
137 1 


133 - 1 
121^3 
110 


98-4 
85-6 

751 


80-8 
66^5 
55^8 


61-5 
4.1-6 
315 


•38 


•385 


3720 


18Jf.-l 
168^0 
162^0 


171-9 
156^0 
140 9 


141-0 
126-1 
111^6 


103-8 
87-9 
75-2 


85-1 
61-1 
546 


65-1 
4:1-0 
29-5 


•50 


•357 


2960 


205 • 
178^0 
151 7 


191-5 
165^0 
140 


156-4 

i3ro 

108-8 


115-3 
89-6 
695 


94'5 
66-2 
476 


73-1 
36-4 
225 


•60" 


•341 


2140 


185-5 
153^0 
123 8 


173-1 

i4r4 

113 7 


Ul-8 

111-3 

87-0 


103-4 
13-0 
52-8 


84-0 
52^1 
343 


68-0 
26-b 
161 

i 



Relative Intensities from Zenith Sky Calculated frvm Mount Wilson Observations in Terms of Mean Attenuation 

Coefficients. 



T(0,0 = ^{f(l+cosH)}eS 
VOL. COXII. — A. 



Ce-'^Q-{C(sec {-!)} + J^h ) * (0, 0) 
3 G 
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Table Yll.— -Total Solar and Sky Radiation Calculated at Mount Wilson Level from Mean 

Attennation Coefficients. 

In order to obtain tlie total solar intensity reaching the earth's surface for different zenith distances of 

the sun, a table of the values of Se~^^^^^ was drawn up for the various values of 0, X, S, and {"given in 
Table VI. From these curves were drawn on a large scale, which, when integrated, gave the values of 



'00 



'00 



E(O^X- Se~^''''^dX, 
Jo 

the units of intensity being those given for S in Table I., and the unit of wave-length being taken as 1//, 
and dX = • 1/^. As a check on the calculations it is interesting to see how far the total intensity can be 
represented by a formula of the type 



ro3 



•00 



E{OdX = ^-csecn g^^x. 



The values of the " apparent " coefficient of attenuation calculated from the integrated area of each of 
the curves of solar intensity corresponding to different values of ( are given below :— 



c 


0°. 


20°. 


40\ 


60°. 


70^ 


80^ 


Mean values C between 
0° and 60°. 





•1023 


•lUO 


•1293 


•1079 


•0953 


•0817 


= -1134,6-^= •893. 



■c _ 



895, 



The value of the apparent transmission of total radiation given for Mount Wilson"^ is e' 

C = •111, in good agreement with the above values. It will be noticed that an exponential formula for 
the total intensity fails for greater zenith distances than 60 °e 

From the intensities from zenith sky given in Table VI., and drawn in Diagrams II.-VIL, the integrals 

T (0, () dX were calculated for various values of {". From these the total sky radiation on a horizontal 

Jo 

plane, H (Q, was calculated from the approximate formula (68), and is plotted against values of ( in 
iagram IX. It will be noticed that the value of the ratio 



D^- 



'00 



H (0 dxj 



roD 



E (0 d\ 



is very nearly constant for all zenith distances and has the mean value '050. This value is in fair 
agreement with the value ^052 (August 18, 1905) but is considerably smaller than the value ^077 
(September 8, 1906, October 19, 1906) determined at Mount Wilson.f It will be noticed under 
Table XII. that the attenuation coefficients for August 18, 1905, are in much better agreement with the 
mean values employed in these calculations than those for October 19, 1906. 

The entries in the last row of the present table give the ratio shy/sun for equal solid angles, the semi- 
diameter of the sun being taken as 16' of arc. These values are those most often measured in observations 
on sky radiation ; the results are shown graphically in Diagram XL 



* ' Annals,' vol, II., p. 96. 

t ' Annals,' vol. II,, Table 35, p. 153. 
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Table YII. 



Units. 



0^ 



20°. 



40°. 



60°. 



70°. 



Total Solar Eadiation per unit area 
normal to sun^s rays reaching 
earth's surface. 

Values of 

Jo 

Too •' 

I Bd\ = 3S21 (S, X) units - 1 • 922 
calories per sq. cm. per minute. 



(S, X) units 
dX == * IfjL. 



3000 • 5 



2944-5 



2782-5 



2673-0 



Calories 
per square 
centimetre 
per minute. 



1-740 



1-706 



1-623 



1-550 



Total Solar Eadiation per unit area 
of horizontal plane reaching earth's 
surface. 

Values of 

GOB C {^ 'E (0 dX 

JO 



(S, X) units 
dX = - 1/x. 



3000 • 5 



2770 



2130 



1337 



Calories 
per square 
centimetre 
per minute. 



1-740 



1-602 



1-242 



775 



Total Eadiation per unit area of hori- 
zontal surface from unit solid angle 
of zenith sky. 

Values of 







T (0, dX. 



(S, X) units 
dk = - 1/x. 



SS - 85 
36-15 
38-05 



S2'2 
34-4 
36-6 



25-5 
27-8 
30-0 



18-3 
20-3 
22-6 



Calories 
per square 
centimetre 
per minute. 



0196 
0209 
0220 



0186 
0199 
0212 



0147 
0161 
0173 



0106 
0117 
0131 



Total Sky Eadiation per unit area of 

horizontal surface. 
Values of 



^^^ f(l+COs2 0Jo 



00 

T{0,OdX. 



(S, X) units 
dX = • IfjL. 



Ul'9 
151-2 
159-3 






US' 

153-1 
163-0 



134' 8 
147-0 
158-6 



2 



122' 

135-9 

151-1 



Calories 
per square 
centimetre 
per minute. 



'0824 
'0S11 
-0925 



0832 
0888 
0945 



'0782 
-0853 
•0920 



0710 
0788 
0876 



Total Eadiation on horizontal surface 
'{cosCE(0 + H(0}^^X. 



Calories 
per square 
centimetre 
per minute. 



1'822 
1-82 



Q 

KJ 



1-833 



1'685 
1-691 
1-697 



1'320 
1-327 
1-334 



'846 
•854 
-863 



■00 

B.{i)dX 



Eatio 



'OD 



cosn '¥j{i)dX 



'047 3 
-0504 
•0531 



0517 
0553 
0588 



0633 
0690 
0744 



0915 
1018 
1131 



2512-0 



1-457 



858 



498 



13'4 
15-9 

18-6 



0077 
0092 
0107 



100 • 6 
119-2 
139-6 



- 0584 
-0692 
-0809 



556 
567 
597 



Eatio 



H(Of/X 
E (0 dX 



' 0473 
-0504 
•0531 



0487 
0520 
0553 



'0485 
•0529 
•0570 



• 0453 
•0508 
•0665 



1171 
1390 
1623 



0400 
0475 
0555 



Eatio (0 



foo 



d\ 



'00 



(0 



'E{i)dX 
= 27r(l-cosl6') = 6-80x10-5 



76 


•6'x 


10- 


-8 


82 
86 


•0 
•2 







74'3xl0-^ 

79-5 

84-5 



62 


'2x 


i6>- 


-b' 


67 


•9 






73 


2 







46'6xlO-^ 

5r7 

575 



.f^' 


3x 


i(9- 


-8 


43' 


1 






50' 


3 







80' 



2018-5 



1-169 



350 



203 



9-8 
11-4 
14-9 



0066 
0086 



79 '4 

92^5 

120 9 



'O46I 
•0537 
•0701 



•257 
•273 



227 
264 
346 



• 0394 
•0458 
•0597 



33' Ox 10-^ 

38^4 

502 



Total Solar and Shy Radiation Calculated at Mount Wilson Level from Mean Attenuation Goefficients. 

3 a 2 
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Table YllL—Eelatim Intensities from Zenith Sky Calculated from Washington Observations in Terms of 

Mean Attenuation Coefficients. 

TMs table is constructed on tlie same plan as Table YI. for Mount Wilson, making use of the attenuation 
constants given in Table II. for WasMngton. The effect of *' dust '' requires us to make use of different 
constants for long- and for short-wave radiation in calculating wave-lengths corresponding to the tabulated 
values of C from the formula G ^= pX~^ + y. The result is a discontinuity in the sky radiation curve in 
the neighbourhood of -GlO/x, shown in the curves of Diagram XII. The mean value of the solution of 
the integral equation is alone given, since the values of the absorption coefficients vary so rapidly from 
day to day that the intensities of sky radiation calculated from them probably differ by more than the 
difference between the extreme solutions, 
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Table VIII. 



c. 


X. 


s. 


0\ 


20^ 


40^ 


60^ 


70°. 


80^ 

r25 


•10 


1^08 fx 


1360 


3-04 


2^85 


2-37 


rai 


l'5A 


•u 


•824 


2470 


16^2 


15-2 


12^5 
29-9 


9^44 
22^15 


7-85 


6^11 


•18 


•725 


3350 


38-5 


36^1 


18^23 


13^4 


•22 


•654 


4240 


69-0 


64^9 


52^8 


38^70 


3r40 


22-1 


,•26 


•625 


4640 


97^5 


91 •I 


74^4 


53^8 


42 • 9 


29^3 


•30 


•556 


5460 


58^5 


54^5 


44^5 


32^1 


25^6 


16^0 


•34 
•38 


•509 


6020 


88^0 
127^1 


83^0 


67^1 


473 


36 • 4 


22^8 


•482 


6230 


109^0 


88-0 


61 •I 


46^5 


28-4 


•50 


•426 


5320 


163^1 


150^9 


119^4 


80^5 


59^4 


32-8 


•60 


•398 


4320 


173^0 


160^0 


126^0 ' 


82^4 


57^3 


30-8 


•70 


•376 


3460 

■ 

2900 


169^8 


156^0 


121^8 


77-5 


53^9 


27^0 


•80 
•90 


•360 
•346 


166^9 


153-1 


118^2 


73^6 


49^7 


.. 22^8 


2240 


146^0 


133-8 


102^6 


62^0 


41-3 1 


19^5 



Relative Intensities from Zenith Sky Calculated from Washington Observations in Terms of Mean 

Attenuation Coefficients. 

T(0.0 = ^{f(l+cosn)}gSrC«-<=G{C(secf-l)}+J^.*(C,0)' 



414 



MR. LOUIS YESSOT KING ON THE SCATTERING AND 



Table IX.— fatal Solar and Shy Radiation Calculated, at Washington Level from Mean Attenuation 

Coefficients. 

This table is drawn up on the same plan as Table VII., making use of the Washington data given in 

Table II. The values of the " apparent " coefficients of attenuation, C, derived from the integrated solar 
radiation curves corresponding to different zenith distances of the sun are given below :— 



C 



C 



o^ 



247 



20°. 



247 



40°. 



•252 



60°. 



'246 



70°. 



•235 



80°. 



210 



Mean value G between 
0° and 60°. 



- -248, e 



-G 



780 



The value of the apparent transmission of total radiation given for Washington* is e~^ 

C = • 240 in substantial agreement with the above mean value between 0° and 60°. 
The ratio 



787, 



H(0^x/r E(C)JX 



is approximately constant for all zenith distances to 70° and has the mean value -078. The value of this 
ratio f or ^ = 80° seems to indicate the existence of a rapidly increasing value of the ratio beyond that 
angle : the analysis by which sky radiation is calculated ceases to hold even approximately for greater 
zenith distances on account of the curvature of the earth. 



* * Annals/ vol. II., p. 96. 
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Table IX. 



Units. 



0\ 



20^ 



40^ 



60\ 



Total Solar Eadiation per unit area 
normal to sun's rays reaching 
earth's surface. 

Values of 

Ceo 

SdX = 3321 (S, X) units = 1 • 922 
calories per sq. cm. per minute. 



(S, X) units 
dX = ' 1/x. 



2597-5 



2553-0 



2389-0 



2032 • 5 



Calories 
per square 
centimetre 
per minute, 



1-502 



1-478 



1 - 382 



1-177 



Total Solar Eadiation per unit area 
of horizontal plane reaching earth's 
surface. 

Values of 



'CD 



cosfl E(O^X. 





(S, X) units 
dX = ' 1/x. 



2597-5 



2400 



1830 



1016 



Calories 
per square 
centimetre 
per minute. 



1-502 



1-390 



1-060 



589 



7o^ 



1670-0 



967 



571 



331 



80\ 



991-0 



574 



172 



100 



Total Eadiation per unit area of hori- 
zontal surface from unit solid angle 
of zenith sky. 

Values of 



T (0, dX. 



(S, X) units 
dX = • Ifjb. 



47-5 



44-1 



34-9 



24-1 



Calories 
per square 
centimetre 
per minute. 



0275 



0255 



0202 



-0140 



Total Sky Eadiation per unit area of 

horizontal surface. 
Values of 



^RiOdX^ 





27r 



Kl+cos^O Jo 



•CD 

T{0,OdX. 



(S, X) units 
dX = ' 1/x. 



199 



197 



184 



161 



Calories 
per square 
centimetre 
per minute. 



115 



114 



1065 



0931 



Total Eadiation on horizontal surface 

D 

{cosCE(0 + H(0}^X. 



Calories 
per square 
centimetre 
per minute 



1-617 



1-504 



1-166 



682 



Eatio 



'03 



H (0 dX 



cost I E(O^X 



0778 



0821 



1005 



1588 



17-4 



0101 



130 



0753 



406 



228 



10-6 



0062 



86 



0497 



150 



506 



Eatio 



E (0 dX 



0778 



0773 



0766 



0793 



1(0, Od\ 



Eatio 0) 



0) 



^yiodx 

27r(l-cosl6') = 6-80x10-5. 



124-3x10-8 



117-6x10-8 



99-4x10-8 



80-7x10-8 



0779 



70-8x10-8 



0865 



72-7x10-8 



Total Solar and Sky Eadiation Calculated at Sea-Level from Mean Attenuation Coefficients at Washington. 
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Table X. — Quality of Skp Badiation at Mount Wilson, October 17, 1906. 

Observations on the quality of sky-radiation have been made at Mount Wilson.* The values given in 
the present table are calculated in terms of the mean coefficients of attenuation at that station. The units 
of intensity in the calculated and observed values aje both arbitrary. The calculated intensities are 
reduced to the same units as the observed intensities by multiplying the former by the mean value of the 
ratios of the intensities corresponding to the different wave-lengths. The results are compared graphically 
in Diagram XIII. The agreement is seen to be fairly satisfactory, in spite of the fact that the msan values 
of the coefficients of attenuation are used in calculating intensities on a particular day. 



'Annals, vol. il., p. 155. 
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Table X. 





Solar 

intensity, 

C = 43°. 


• 422/x. 
3430 


• 457/x. 
4700 


1 

•491/x. 

4960 


•556//.. 
4630 


•614/x. 
4080 


• 660/x. 
3600 


(i-) 
Quality of Sky Radiation. 

Comparison of 

Calculated Values with 

Mount Wilson 

Observations, 

October 17, 1906. 

^ = 23°. 

xp = 79°. 


Zenith 

sky, 
C = 43°, 


lU 

116 
109 


102 
97 
92 


88 
80 
78 


U 

43 
42 


26 
26 


18 
18 


Skyjsun, 


361 X 10"^ 

339 

318 


218 X 10-^ 

206 

196 


167 X 10"^ 

161 

157 


95-0x10-4 

92-5 

90-5 


68-5x10-4 

63-5 

63-5 


50 X 10-4 

50 

50 


Observed 

sky/sun. 


655 


521 


294 


188 


106 


100 


Calculated 

sky/sun. 


72% 
678 
636 


^86 
412 
392 


884 
321 
314 


iP(9 

185 
181 


127 - 

127 

127 


100 
100 
100 


(ii.) 

Quality of Sky Radiation. 

Comparison of 

Calculated Values with 

Mount Wilson 

Observations. 

C = 58°. 
c/) = 17°. 
ij^ - 146°. 

(iii.) 

Quality of Sky Radiation. 

Comparison of 

Calculated Values of 

Zenith Sky with Mean 

Absolute Intensity of 

Sky Eadiation. 

Mean, C = 50°. 


Solar 

intensity, 

t - 58°. 

Zenith 

sky. 

t = 58°. 


3050 


4210 


4570 


4360 


3850 


3420 


99-0 
91-0 
82-5 


81-0 
77-5 
73-5 


65 
62 
60 


87 
36 
35 


20-5 
20-5 


18-5 
13-5 
13-5 


Sky/sun, 
C = 58°. 


S25 X 10-'^ 

295 

271 


198 X 10-^ 

184 

175 


14^2 X 10-4 

136 

131 


85-0x10-4 

82-5 

800 


58-8x10-4 

53-3 

53-3 


89-5x10-4 

39-5 
39-5 


Observed 

sky/sun. 

Calculated 

sky/sun. 

Zenith 

sky. 

Observed 

mean 

sky. 


574 


425 


317 


191 


124 


104 


750 
680 
625 

112-6 

104-5 

970 


U5 
424 
403 

92-0 
88-0 
83-6 


827 
314 
302 

78-5 
71-5 
690 


196 
190 
185 

40-5. 

40-0 

39-5 


123 
123 

22-2 
22-2 
22-2 


91 
91 
91 


15-7 
15-7 
15-7 


1194 


986 


701 


395 


231 


174 


1 
Calculated 
zenith 
sky. 


1192 
1109 
1030 


975 
933 
885 


780 
758 
731 


429 
424 
418 


285 
235 
235 


166 
166 
166 
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Table XL— Total Shy Badiation at Motmt fViho% October 19^ 1906. 

Observations on the total intensity of sky-radiation are given for Mount Wilson* under the date 
October 19, 1906. Attenuation coefficients for this day are not given, m that the observations are 
compared ^vith values calculated from the mean coefficients at Mount Wilson. If w is the solid angle 
subtended by the sun, the value of 



■00 

CO 



is calculated from the approximate formula (^^J^ 



m{T{cj,,Od^ = <oi±-5^seC(/>|'^T(0,0^X. 

I + COS"^ ( Jo 



Making use of the values of 

* TOO 

T (0, dX / E (0 dX 



'05 



H) 



from Table VII., and interpolating for intermediate values of {" from the curves given in Diagram XL, the 
ratio 

'oo I Too 

was calculated. 

The angle between sun and sky was calculated from the formula 

cos = Qo&cj) cos {■+ sin ^ sin f cos f, 

the azimuth \p being measured from a vertical plane through the sun in the direction N.-E.~S.~W. 

The comparison of the calculated ratios with those obser\red is fairly satisfactory ; the greatest 
discrepancies occur when is small and when {" is large. The probable reasons for this lack of agreement 
are given at the end of Section 8. 



* ' Annals/ vol, II., Table 82, p. 151 
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Table XI. 



i 

! 
i 


c/>. 


^. 


e. 


1 + cos'^ 

1 + C0S2 ^' 


sec c^. 


Zenith sky ^ -^^.g 
Sun 


f^xlO-B. 
8UE 


Observed 
Sky/Sun x 10"^. 


55-8 


79-2 


_^ 43-5 


d'l 


1*127 


5-34 


60-2 
55*0 
49-0 


861 
330 
294 


295 


55-5 


67-4 


- 41-2 


37-9 


1-228 


2-60 


60-6 
55*1 
50-0 


192 
176 
160 


219 


54-0 


54 -e^ 


- 35*5 


28-7 


1-316 


1-72 


61-8 
56-3 
50-5 


HO 

127 
114 


165 


53-6 


44-6 


- 21*3 


18-6 


1 - 402 


1-40 


6fd ' 
57-0 
510 


112 
100 


190 


50-4 


40-5 


- 6-7 


11-5 


1*392 


1-31 


6^-6 
59-3 
53-2 


118 

108 

97 

101 
93 
84 


285 


49*8 


31*1 


12*9 


20*6 


1*327 


1-17 
1 • 55 
1^57 


65-0 
60-0 
54- 


350 


49*2 


■ 

49*7 


-149*6 
™ 113-3 


94-4 
79-1 


•696 


65-5 
60 • 4 
54-5 

6o'8 
60-6 
54-8 


70-7 
65^2 
58-8 

7S'0 
67-2 
60-8 


98 


48-9 


50-4 


-707 


108 


47^4 


79-2 


- 63-4 


63 - 2 


•826 


5-34 


67-0 
61-8 
560 


%96 
273 
247 


222 


47-0 
46-5 


77-9 
29*2 


- 84-5 


77'- 7 


•714 


4-77 


67'% 
62-0 
56-2 

61-8 
62-5 
570 


m9 

211 
191 

76 ' 1 
70-2 
64- 


171 


- 79-6 


48-5 


•977 


1*15 


98 


46-2 


27-8 


144-8 


70-3 


•752 


1-13 


68-0 
63-0 
57-1 


S7'7 
53*5 
48-5 


68 



'Total Sky Radiation at Mount Wilson, October 19, J 906. 
Comparison of ratio shjjsun for equal solid angles with results of calculation from the formula 



(a 



{i,i)dxl 



°E (0 d\. 



O H. 4i 
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Table XlL—^-Ammge Intsnsiiies oj Sky Madiation, Mount Wiisork 

Observations are given at Mount Wilson* for the average total intensity of sky radiation taken for 
different azimntlis. The theoretical formula corresponding to average' sky intensity is by {%%\ 



1 



9, 



w 



'2Tr 



CO 







# Ti<l>,OdX 



27r 



— ..^ — ^ dip 
1 + cos^ I 



sec </j 



'00 

T (0, dX, 



The term in square brackets denoted by 



1 +(iOB^0^ 
1 H-cos^ {■ 



av. 



1 - J sin2 cj> . 



3 cos^ t - 1 
1 + cos^ {" 



In calculating average shjjsun for equal solid angles, the value of 



CO 



'00 



'00 



^ T(0,OclA/l E{OtlA 



is derived by interpolation from Diagram XL 

In order to compare the attenuation coelhcients on the dates of the above sky observations, these values 
taken from, the Mount Wilson observationsf are compared in a table with the mean values employed in 
the calculation. It will be noticed that the absorption constants on August 18, 1905, are in much better 
agreement with the mean values than those of September 8, 1906. The result is that as we should expect 
the sky observations on the former date agree more closely with the calculated values than those on the 
latter. 



*' 'Annals/ vol. II., Table 32, x>. 151. 
t * Annals/ yoL II., pp. 96-97, 
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Table XII. 
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Sky Radiation, Mount Wilson, September 8, and October 19, 1906. 
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Skj Eadiation, Mount Wilson, August 18, 1905. 
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AUenxiation Coefficients on Dates of Sky Observations, 
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Table XIIL- — Polarization of Sky Radiation Calculated from Mount Wilson and Washington Observations, 

Numerical results on the polarization of sky radiation are calculated from the approximate formulae 
given in Section 6. In order to save space numerical values corresponding to tlie formulae given at the 
foot of the table are entered in the same order under each value of X and {". 

The values of 

P(0,f) - Ce-%{C(secf-»l)} 

are obtained from Tables VI. and VIII. 

The ratios of the intensities of zenith sky polarized in the principal plane and in a plane at right angles 

to it are calculated in the form 

Ti (0, i) „ l-fQ(O.D/P(0,D . 
T2{0,0 cos2C+Q(0,0/P(0,0 

The results are shown graphically in Diagram XIV. 

The ratios Q (0, {")/? (0, {") given in the present table compare the contribution of self-illumination to the 
sky intensity of dilFerent wave-lengths with the effect of direct sunlight on sky radiation. 
The ratios 

Q(o. D/P(Q.D 

1 + Q (0, o/p (0, 

compare the residiml intensity of sky radiation over the great circle of most complete polarization 90° from 
the sun to the intensity polarized in a plane at right angles. It will be noticed that this ratio remains 
nearly the same for all zenith distances, and is very small for long wave-lengths when the effect of self- 
illumination is small. 
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Polarization of Sky Radiation Calculated from Mount Wilson and Washington Observations, 
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Diagram L — VariaMon of Coefficients of Afmospheric Attenuation ivith JFcwe- Length, 

The curves of Diagram I. are drawn from the data given in Table I. The coefficients of attenuation Co; 
at the various stations are plotted against the inverse fourth power of the wave-lengths, the wave-lengths A 
being expressed in microns. The observations give rise to a law of the form 

The stations to which the various lines refer are as follows :-~™ 

Curve I. Mount Whitney. Curve III. Potsdam. 

Curve II. Mount Wilson* Curve IV. Washington. 

The conclusions derived from these observations are described in Section 7, and the constants of 
atmospheric scattering and absorption derived from similar curves plotted on a large scale are given in 
xaDie J-J.. 

The discontinuity in the curves in the neighbourhood of '610//- at low-level stations is clearly marked, 
and the observations show that the curves may be represented by broken straight lines. We cannot 
suppose, however, that the discontinuity at X = -BlOfo will in reality be as sharp as the intersection of 
two straight lines would require. The precise cause of this well-marked point of discontinuity furnishes 
a point which must be left to further investigation. 
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Diagram I. Variation of coefficients of atmospheric attenuation with wave-length. 
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Diagrams II.-VII. — Relative Intensities from Zenith Sky Calculated for Different Wave-Lengths at Mount 

Wilson from Mean Coefficients of Attenuation, 

The curves given on Diagrams II.-VII. are calculated for various wave-lengths from formula (64) in 
terms of the mean coefficients of attenuation for Mount Wilson given in Table II. The numerical values 
from which these curves are drawn are given in Table YI. The three outer curves of each set represent 
the extreme and mean solutions of the integral equation (64)5 

T (0, = (47r)-i [f (1 + cos^ 0] S [P (0, + Q (0, 0] • c/C, 
where 

P(0,0 == Cg"^G[C(secC-l)] 
and 

Q (0, = i g C^\ '^ (C, 0). 

The inmost curve of each set gives the intensity from zenith sky neglecting self-illumination, i.e., it 

represents the value of 

(4^r[i(l+cosH)]SP (0,0.^/0. 

In this way the contribution of self -illumination to the intensity of sky radiation is made clear. We 
notice the divergence of the extreme solutions for large values of the sun's zenith distance. 

The curves just described resemble the curves obtained by NiCOLS* in his observations on the intensities 
in the spectrum of zenith sky compared with corresponding intensities in the spectrum of an acetylene 
flame. The results given by NiCOLS were taken at various times of the day (July 18, 1907) at Sterzing 
in the Tyrol (lat. 46° M N., long. IT 25' E., elevation 3110 feet above sea-level). Numerous observations 
have been taken by Crova,! Zettwuch|, Majorana,§ and others.|| The results are not strictly 
comparable with theory since the attenuation coefficients at the time of observation are not given. The 
forthcoming observations of the Smithsonian Astrophysical Observatory on sky radiation at Mount 
WhitneylT (August, 1910) will enable an accurate comparison to be made with the results calculated from 
the attenuation coefficients measured on the same day. 



* NiCOLS, E, L., '' Theories of the Colour of the Sky," ' Physical Eeyiew,' XXTI., June, 1908, p. 497. An extensive 
bibUography on sky obseryations is given in this paper. 

t Oeova, 0. E., ' Comptes Rendns,' CIX., p. 493 j OXII., p. 1178 ; also ' Annales de Ohimie et de Physique,' (6), XX., p. 480. 

X Zbttwuch, G-., " Ricerche snl Bleu del Oielo," translated in the ' Phil. Mag.,' (6), lY., August, 1902, p. 199, 

§ Majorana, Q., " On the Relative Luminous Intensities of Sun and Sky," ' Phil. Mag.,' (6), I., May, 1901, p. 555. 

[II The writer is indebted to Dr. Otto Klotz, of the Dominion Observatory, Ottawa, for calling his attention to the work of 
ExNER along these lines. Exnbe, ' Sitzungsbericht d. K. Akad. d. Wissen., Wien, M. N. Klasse,' 1909, vol. 118, IIa. A 
summary of Exnbe's work and of the observations of Wiisnkr and Schramm, W., is given by Abbot, 0. G-., in his recent 
book, * The Sun,' p. 299 (Appleton and Co., 1911). A problem of somewhat the same nature as the present one is considered 
by G-OLB, E., "The Isothermal Layer of the Atmosphere and Atmospheric Radiation," 'Roy. Soc. Proc.,' A, vol, 82, 
im^.—Note added December 31, 1912.] 

^ Abbot, C. Q-., '' Report on the Astrophysical Observatory," ' Annual Report of the Smithsonian Institution,' 1911, 

p. 65. 
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Eelative intensities from zenith sky calculated from Mount 
Wilson observations. ^ = 0°. 
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Wilson observations. 
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Diagram IV. Eelative intensities from zenith sky calculated from Mount 

Wilson observations. {" = 40°. 
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Diagram V. 
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Eelative intensities from zenith sky calculated from Mount 
Wilson observations. ^ = 60°. 



Diagrams il.-VII. 
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Diagram VI. Relative intensities from zenith sky calculated from Mount 

Wilson observations. {" = 70°. 
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Diagram VII. Eelative intensities from zenith sky calculated from Mount 

Wilson observations. ( = 80°. 
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Diagram VIII. — Comparison of Formulm far the Relative Intensities of Sky Radiation Calculated from Mount 

Wilson Observations for {" = 60°. 

It is of some interest to compare with the formula for sky intensity derived in the present paper the 
various formulae which have been employed by various writers in the interpretation of sky observations. 
The results are shown graphically in Diagram VIII. for a zenith distance of the sun of 60°, and refer to 
intensities of zenith sky. 



Curve I. is given by the formula (63) 



(47r)--l Sf (1 + C0S2 C, 



and is the same as that given by Kelvin^ and employed by PEURlNt in an estimate of the 

number of molecules per unit volume of a gas. 

Ctirve 11. is given by the formula 

(47r)-iS{|(l+cos2 0}c, 

and is seen to give a considerably smaller value of sky intensity in the red, agreeing closely with 

the value T (0, f) for wave-lengths greater than -45/*. 
Curve III. is given by the mean solution (64) for T(0, C) which takes into account the effect of self- 

illumination. 
Curve IF. gives the intensity of zenith sky neglecting self -illumination but taking into consideration the 

attenuation of both the incident and scattered radiation. 



* Kelvin, 'Baltimore Lectures' (1904), p. 313. 

t Perrif, J., * Annales de Ohimie et de Physique,' S'"^ Serie, September, 1909. 
Molecular Reality' (F. Soddy, 1910), published by Taylor and Francis, London. 



See p. 79 of ' Brownian MoTement and 
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Diagram YIII. Comparison of formulae for intensities from zenith sky, Mount Wilson, ( = 60". 
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DiAGBAMS IX. AND \.— Total Radiation Calculated from Mount Wilson and Washington Observations. 

The curves drawn on Diagrams IX. and X. are drawn from the data calculated in Tables VII. and IX. 
for the values of total solar and sky radiation at Mount Wilson and Washington respectively. The 
various curves denoted by I., II., III., IV., refer to the following quantities: — 

Curve I. gives the total intensity of solar radiation in calories per square centimetre per minute normal 
to the sun's rays ; i.e., the curve represents the variation of 

E (0 d\ with i. 

Curve //.gives the total intensity of solar radiation in calories per square centimetre per minute on a 
horizontal plane; i.e., the curve represents the variation of 



■00 



cosq E(0^^witht. 

Curve III. represents the total intensity in calories per square centimetre per minute of the radiation 
from the entire skj^ on a horizontal plane; i.e., the curve represents the variation of 



[" H (0 d\ with t. 



Curve IV. represents the total radiation from sun and sky on a horizontal plane ; i.e., the curve 
represents the variation of 

'""' {cos CE (0 + H(0}rfX with C. 







The calculations are only given as far as f == 80°. Curves I. and II. fall to zero at {" == 90°. Curve III. 
is arbitrarily carried on beyond this point (in dotted lines) to f = 100°, and represents to some extent the 
intensity of scattered radiation due to the " twilight " following the sunset or preceding the sunrise. The 
effect of the earth's curvature makes it practically impossible to calculate the scattered radiation when the 
sun is below the horizon. The heat reaching the earth before sunrise shows itself in continuous air- 
temperature records in the fact that the minimum in the diurnal inequality of temperature occurs about 
one-half hour before the sun actually rises. 

The curve just described are intended for use in meteorological problems. 
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Diagram IX, Total solar and sky radiation calculated from Mount Wilson observations. 
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Diagram X. Total solar and sky radiation calculated from Washington observations. 
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Diagram XL — Relative Intensities from Sky and Sun Calculated, from Mount Wilson and 

Washington Observations. 

The curves given in Diagram XL represent the ratio zenith shyjsun for equal solid angles, the total 
radiation covering all wave-lengths being taken in both cases. The curve represents the formula 



(0 

'o 



T (0, f^l r E (0 dk, 



where w is the solid angle subtended by the sun, ^.^., (o ^ 27r(l - cos 16') = 6*80 x 10 ^. 

Curves I. refer to Mount Wilson and include the extreme and m.ean solutions in the expression for 
T (0, 0- 

Curve II. refers to Washington and represents the mean solution of the integral equation only. 

Since the ratio just described is that most easily obtained in measurements of sky radiation, the curves 
just given are appropriate for a comparison of theory with observation. 



Diagram XIL — Relative Intensities from Zenith Sky, Calculated for Different 'Wave-Lengths from 

Washington Observations. 

The curves given in Diagram XII. are drawn on the same plan as those described in Diagrams IL-YII. 
for Mount Wilson. In the present instance the numerical values from which the curves are drawn are 
given in Table VIII. Intensities corresponding to the mean solution of the integral equation are alone 
given, since the presence of " dust " at lowdevel stations is an extremely variable factor. This factor gives 
rise to a sharp bend in the attenuation curve for Washington shown in Diagram L in the neighbourhood 
of -eiO/x. The effect on the intensity of zenith sky is very sharply marked and gives rise to a finite 
discontinuity in the intensity curve corresponding to that wave-length. In reality the bend of the 
attenuation curve in Diagram I. cannot be as sharp as the intersection of two straight lines, so that the 
discontinuity in the sky intensity curve in this neighbourhood is probably represented by a sharp peak 
(represented in the drawing by the dotted portions of the curves) which arises from the effect of '' dust '^ in 
scattering the long-wave radiation. The prominence of this peak in the intensity curve must be an 
extremely variable factor, and may possibly account for discrepancies in observations on sky-radiation ; e.g., 
NlCOLS"^ mentions as a result of his measurements that sky intensities show far greater relative intensities 
in the longer wave-lengths than one would expect from Rayleigh's theory. The effect of *' dust " is an 
illustration of the effect of several groups of molecular complexes or small particles on radiation travelling 
through a medium containing them ; each group which is able to produce a sharp bend in the attenuation- 
curve at any wave-length is also able to produce a corresponding maximum in the intensity-curve of the 
scattered radiation, the whole effect resembling a broad emission band in its spectrum, f 

* NicoLS, loc. clt., p. 502. 

[t Important experiments describing the transition from lateral scattering to selective reflection are described hy 
Wood, R. W., under the name "resonance radiation." These are described in the 'Clark University Lectures,' Clark 
University, 1912, p. 101. — ^ote added December 31, 1912.] 
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Diagram XL — Eatio zenith sky /sun for equal solid angles calculated from Mount Wilson and 

Washington observations. 
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Diagram XII. Eelative intensities from zenith, sky, calculated from Washington observations. 
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Diagram XIII. — Quality of Shy Radiation Calmlated at Mount Wilson and Compared with Observations. 

The curves (i), (ii), (iii) of Diagram XIII. are plotted from the results of Table X., and show graphically 
how observed values of sky radiation agree with values calculated from mean attenuation coefficients at 
Mount Wilson. The curves corresponding to the extreme and mean values of the solution of the integral 
equation are retained throughout. Observed points are denoted by small circles, and lie fairly close to the 
calculated curves. 



Diagram XIV. — Polarization of Zenith Sky Calculated from Washington and Mount Wilson Observations. 

The curves given in Diagram XI Y. are drawn from data calculated in Table XIIL, and represent the 
ratio of the component intensities of zenith sky polarized in the principal plane to the component polarized 
at right angles to that plane. This ratio is represented roughly by the formula 

Ti(o,0 __ 1 + Q (0, 0/P (0, C) 



T2(0,0 cos2f+Q(0,0/P(0,0 

Curves are given for various zenith distances of the sun, making use of the constants for Mount Wilson 
and Washington. In the latter case the effect of " dust " gives rise to a discontinuity in the neighbourhood 
of ' QIO fjb, represented in an actual case by a peak in the curves (represented by dotted lines). The 
characteristic just mentioned is actually observed in curves of polarization of zenith sky obtained by 
NlCOLS."^ In comparatively " dust-free " air as at Mount Wilson this peak is wanting, and the curve 
represents a type obtained by NlCOLSf which rises with increasing wave-length. Since the forms of the 
curves just described depend essentially on the effects of self -illumination, these will be extremely 
sensitive to the presence of " dust," and will therefore vary greatly with time and place. Eeflection of 
solar radiation from the earth's surface (e.g.^ a snow-covered landscape) will affect in a marked manner 
the polarization of sky radiation independently of effects of self-illumination, and may therefore account 
for the great variety of results obtained by, observation on this point. 



* NiooLS, loc. cit., p. 508, fig. 10, Curves (b) and (&'). 
t NiooiiS, loc. cit.f fig. 10, Curye (c). 
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Diagram XIII. Quality of sky radiation, calculated from Mount Wilson observations, October 17, 1906. 




Mount Wilson 



Washington 



Diagram XIV. Polarization of zenith sky calculated from Mount Wilson and Washington observations. 
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